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A METHOD OF ASYMPTOTIC EXPANSIONS FOR SINGULAR 


PERTURBATION PROBLEMS WITH APPLICATION 

IN VISCOUS FLOW 
By E. Dale Martin 

Ames Research Center 
SUMMARY 


A method of asymptotic expansions for solving certain classes of 
singular perturbation problems is presented. The method is a new approach to 
constructing inner and outer expansions. A classification of singular pertur- 
bation problems is given to clarify when this method is applicable and to aid 
in determination of an appropriate matching rule to use. Simple examples are 
given to illustrate use of the method. 

For an important class of boundary-layer-type problems, this approach 
employs a stronger form of the matching principle than used previously. All 
displacement effects are retained explicitly in this matching process, and 
higher order displacements are identified and calculated. The concepts of com- 
plex displacement and multiple displacements are introduced. The development 
of the method leads to clarification of conditions for algebraic or exponen- 
tial approach of inner to outer expansion solutions. 

The expansion method developed here is applied to a model of a familiar 
typical physical problem in viscous -flow theory: hypersonic flow over a 

sphere. This application is viewed as an investigation in viscous flow of the 
role of displacement In matching the boundary-layer asymptotic expansion to 
the outer expansion. Higher order boundary -layer displacements are calculated 
explicitly. Results for skin friction, shock standoff distance, etc., are 
computed and compared with results from a numerical solution of the same model 
problem to demonstrate the validity of the expansion method. 


INTRODUCTION 


The main purpose of this paper is to outline a method for solving certain 
classes of singular perturbation problems including, in particular, an Impor- 
tant class of boundary-layer-type problems. The proposed method is a new 
approach to defining and constructing inner and outer expansions. This 
approach makes possible the use of a stronger form of the matching principle 
than that used previously. It may be especially useful in problems where 
certain displacement effects are important. 

Singular perturbation problems occur in virtually every branch of 
mathematical physics. A singular perturbation problem Is encountered in seek- 
ing a mathematical solution for a limiting value (e.g., zero) of a parameter. 



or in seeking a small perturbation solution, and when the solution so obtained 
is not uniformly valid throughout the domain of the independent variables. 

(For a comprehensive discussion and treatment, see ref. 1.) For example, the 
formal solution may become infinite at some locus of the independent (e.g., 
space) variables. This occurs near a rounded leading edge in the small per- 
turbation solution of thin-airfoil theory (see ref. 2). In other problems 
the nonuniformity may be evidenced by a discontinuity of the solution within 
the region of interest (e.g., the discontinuous pressure across a shock wave 
in inviscid-f low theory) or by failure of the solution to satisfy a boundary 
condition (i.e., by a discontinuity at the boundary). The prime example of 
this type is in the calculation of fluid flow at high Reynolds number over a 
surface, in which the flow is essentially inviscid except for a thin boundary 
layer at the surface. The inviscid-f low solution fails to satisfy the no -slip 
condition at the boundary surface. The boundary -layer theory of Prandtl 
(ref. 3) was developed to cope with the mathematical nonuniformity in the flow 
calculation. Many other singular perturbation problems of this type occur in 
all branches of fluid mechanics (e.g., see ref. l) . The discontinuity or 
failure to satisfy a boundary condition is usually, but not always, due to the 
order of the differential equation, or set of equations, being reduced in the 
process of letting the small parameter become zero. The discontinuity is 
actually an asymptotic representation of a quick transition (rapid variation) 
of one or more dependent variables over a small range of the independent vari- 
able (ref. 4). In the twenty -eighth Josiah Willard Gibbs lecture in 195^v 
Professor Friedrichs discussed the asymptotic phenomena associated with "quick 
transition regions" in many branches of mathematical physics. Still another 
type of singular perturbation problem has a prime example in the Stokes 1 
paradox of viscous flow, in which the approximate solution valid near a body 
for very small Reynolds number cannot approach the uniform-stream condition at 
large distances. 

Rather than seek a simple limit (asymptotic solution) or simple perturba- 
tion solution for the function of interest, one may attempt to construct a 
perturbation expansion (power series) of the solution in terms of finite 
values of the small parameter. Such an expansion would yield higher approxi- 
mations to the solution (if uniformly valid), which are needed for sufficient 
accuracy if the perturbation parameter is not very small. The small perturba- 
tion theory originated by Sir Isaac Newton has been highly developed by many 
others (see discussion by W. R. Sears, ref. 5). Extension of this theory to 
the asymptotic expansion (power series expansion in the small parameter) was 
devised by Poincar^ (see ref. 6). If any of the terms of the expansion are 
not uniformly valid, ways must be found to determine an asymptotic limit func- 
tion that is uniformly valid and to construct a uniformly valid asymptotic 
expansion extending the solution to higher order approximations for sufficient 
accuracy when the parameter is not very small. 

A number of methods have been developed and used in finding uniformly- 
valid asymptotic expansions (refs. 1 and 7 describe various methods). Each 
method has certain attributes and certain types of problems wherein it is most 
advantageously used. The method of strained coordinates (also known as 
Lighthill’s technique and as the PLK method; see ref. 8) was developed and 
described by Lighthill (refs. 9 and 10) . This method is most useful in treat- 
ing problems of the first type described above (having an infinity in the for- 
mal perturbation solution) and in singular characteristic problems (see 
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ref. 11). In 1953, Kuo (ref. 12) provided an extension of Lighthill*s 
technique (see discussion by Tsien, ref. 8) to a problem in first-order 
boundary -layer theory: viscous flow over a finite flat plate. Lighthill’s 

technique was used to deal with the nonuniformity at the leading edge, whereas 
the nonuniformity associated with the entire boundary layer (caused by the 
essential singularity in the complex plane of the reciprocal of Reynolds num- 
ber) was dealt with by a simple stretching of variables, following Prandtl. 
(Kuo*s treatment of that problem has been criticized by Van Dyke (refs. 1 and 
13 ) on the grounds that it omits a concentrated force at the leading edge and 
that a constant in the solution was obtained by an imprecise asymptotic 
evaluation . ) 

Methods for solving singular perturbation problems of the boundary-layer 
type are extensions of Prandtl* s boundary-layer theory, which incorporates a 
special kind of limiting process that gives a uniformly-valid solution near 
the boundary where the no -slip boundary condition is lost in the inviscid-f low 
theory. In 1904, Prandtl (ref. 3 ) recognized that the outer inviscid-f low 
solution could be joined to an inner solution for the nonuniform region, or 
boundary layer, found by properly magnifying the physical variables and by 
considering the relative order of magnitude of the various terms in the Navier- 
Stokes equations in the limit of infinite Reynolds number. Since then, 
,r boundary-layer methods” have been applied to many other problems (see refs. 4, 
8, 14, and 15; and the literature referred to in those papers) . In 193 5 j 
Prandtl (ref. 1 6 ) suggested how to improve the flat -plate boundary -layer solu- 
tion of Blasius (ref. 17) by iteration. This suggestion led to the use of an 
asymptotic expansion by Alden (ref. l8) with the first-order boundary-layer 
solution of Blasius as the first term. The history of that problem and correc- 
tion and further extension of the solution have been treated in detail by 
Goldstein (ref. 19, chap. 8); the solution was further developed and discussed 
by Murray (ref. 20). In all boundary -layer -type problems, the general ideas 
of Prandtl* s boundary-layer theory provide the first term in an asymptotic 
expansion of the solution to the full governing equations. 

A significant advancement in the use of asymptotic expansions in singular 
perturbation problems was made by the introduction and development of a tech- 
nique by Kaplun, Lagerstrom, and Cole (refs. 21-25). This technique, usually 
referred to as the method of ”inner and outer expansions” (but more recently 
as "matched asymptotic expansions"), employs two distinct asymptotic series 
expansions of the solution: an outer expansion valid away from the nonuni- 

formity and an inner expansion valid in the "inner region" of the nonuniform- 
ity. Limiting processes are defined according to which the inner and outer 
expansions must match in their "overlap region of common validity. " The devel- 
opment of this procedure was closely related to some previous work by Latta 
(ref. 26). The method has been further developed and procedures established 
(especially for boundary-layer -type problems) by Van Dyke (refs. 13 and 27 
to 31) • 

A method based on the ideas of inner and outer expansions is developed 
in the present paper. This method, which employs "displacement variables," a 
stronger form of the matching principle, and some new rules for obtaining 
appropriate forms for the expansions, may be advantageous in some problems 
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involving special displacement effects. The particular definitions of inner 
and outer functions used here and the use of displacement variables represent 
a different point of view from the conventional inner and outer expansion 
method. A result of this different point of view is that a more precise 
matching principle is found to apply in certain problems, which is useful in 
determining and understanding displacement effects. Development of the method 
includes introduction of the concepts of complex displacement and multiple 
displacements. Development of the new matching rule also leads to knowledge 
of conditions for which the inner expansion can approach a form of the outer 
expansion term by term exponentially and of the reasons for either exponential 
or algebraic termwise approach of inner solutions to outer solutions. The 
method is formulated in terms of power-series expansions of some form of the 
small parameter, but terms of higher order than those considered in a calcula- 
tion may be of some order other than a power of the small parameter. 

After the method is outlined and discussed in detail, including several 
simple examples to illustrate procedures, it is applied to a simplified model 
of the familiar problem of hypersonic viscous flow over a sphere at high 
Reynolds number. That problem has been studied by a number of investigators, 
but the main purposes in this application will be to illustrate the precise 
determination of displacement effects, to identify and calculate the higher 
order displacements themselves, and possibly to generate an increased under- 
standing, by the point of view of this method, of the role played by displace- 
ment in the relationship between inner and outer solutions. 


CHARACTERISTICS OF SEVERAL CLASSES OF SINGULAR 
PERTURBATION PROBLEMS 


This section discusses several different classes of singular perturbation 
problems and their distinctions. The main purpose is to make clear the place 
of the method discussed below in singular perturbation theory. Since differ- 
ent types of nonuniformities have different characteristics and must be han- 
dled somewhat differently, it is convenient to classify problems according to 
the type of nonuniformity. This classification is necessary to the discussion 
of different types of inner and outer solutions, of the ways the inner solu- 
tion can approach the outer solution, and, consequently, of the matching prin- 
ciples relating inner and outer solutions. Some important aspects of this 
classification, and of the examples given, are developed and discussed more 
fully in later sections. 

Purely for purposes of reference in this study, different types of 
singular perturbation problems will arbitrarily be denoted by class 1, class 2, 
etc. There are classes other than those discussed here and many problems 
overlap the classes, that is, may have the characteristics of two classes at 
the same time . 

The problems used here to illustrate different classes will involve only 
ordinary differential equations, and will be very simple. However, the 
classes of problems these examples represent, and the asymptotic phenomena 
they exemplify, are extremely important. The same classes include problems 
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■with more than one independent variable and with both linear and nonlinear 
equations. The primary use of* expansion methods is in finding analytical 
solutions to those more difficult problems. Understanding of the simpler 
problems can lead ultimately to the successful solution of the more difficult 
problems with at least some of the same essential features. 

In the following, a dependent variable f may be a function of several 
independent variables, including a small parameter e: 

f = f(x,y,z,e) 

The parameter e may be thought of as a complex variable and f as a complex 
function of e, even though e and f will ultimately take on only real values 
in the solution. The theory of complex variables and, in particular, of ana- 
lytic functions can then be used in discussing the behavior of f as e 
varies, for example, as e approaches zero (cf. ref. 1). Since the examples 
that follow will contain only derivatives with respect to y, they will be 
written as ordinary differential equations for f(y) with d( )/dy denoted 

"by ( ) ' • 

The various classes, to be illustrated below, will be distinguished by 
properties I, II, and III, described by the following statements (where either 
l(a) or l(b) describes property I, etc.): 

I. Order of the differential equation 

(a) is not reduced upon letting e -*■ 0. 

(t) is reduced upon letting e -* 0. 

II. A discontinuity or loss of a boundary condition in the formal 

asymptotic solution obtained for e 0 

(a) does not occur, but the terms of the asymptotic expansion are 
singular (infinite) at some y (e.g., y = 0) ; that expansion 
solution is then designated as an "outer expansion." 

(b) occurs because of a singularity in that expansion solution (then 
designated "outer expansion") at (or near, as e 0) the loca- 
tion where the lost condition was to have been applied. (The 
singularity may not be apparent in lower order terms of the 
expansion. ) 

(c) occurs, but the terms of the expansion solution (then designated 
"outer expansion," to be more completely defined later) are 
otherwise regular (i.e., are sectionally continuous and have 
sectionally continuous, or one-sided, derivatives) at the loca- 
tion where the condition is lost (y = 0) . 

(d) occurs and the limit either of the function or of its deriva- 
tives does not exist as e 0, so there is no "outer solution." 

III. An essential singularity (of the function f) with respect to e 
on the plane e = 0 in f, y, e space 
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(a) does not exist ; an appropriate inner solution is equivalent to 
f, and its asymptotic expansion for small e (the inner expan- 
sion, a valid representation for f near y = 0; to be defined) 
approaches algebraically to the outer expansion term by term. 

(b) exists ; an appropriate inner solution function valid near y = 0 
(to be defined) approaches exponentially to the outer solution 
function as y 00 ; but the asymptotic expansion of the inner 
solution approaches algebraically to the outer expansion term 
by term. 

(c) exists ; and an inner solution (a valid representation for f 
near y = 0; to be defined) and its asymptotic expansion for 
small e approach exponentially to a form of the outer solution 
term by term. 

(d) exists ; but there is no "inner region" (boundary layer , or quick 
transition region) to which the nonuniformity is confined, and 
hence no possible matching of an inner solution with the outer 
solution. 


For convenience, the following table lists the properties of the six 
classes to be discussed. 


Class Properties 


1 

la. 

Ila, 

Ilia 

2 

la. 

lib. 

Ilia 

3 

lb. 

lib. 

Ilia 

4 

lb. 

lib. 

Illb 

5 

lb. 

lie. 

IIIc 

6 

lb. 

lid. 

Illd 


Class 1 will denote the general class of problems for f 
statements la, Ila, and Ilia apply. To illustrate this class, 
example problem 1: 


(y + ef)f r + f = 1 
f(l) = 2 


in which 
consider 


(i) 


This is a special case of a problem used by Lighthill (ref. 10) in illustrat- 
ing the method of strained coordinates and is identical to one used by 
Van Dyke (ref. l) . To study the characteristics of this problem, we note that 
the exact solution is 


f = 


i ± y 


+ 2 ( + 4 


1/2 


_ Z 

€ 


which may be expanded, using the binomial theorem, for small € as: 


f S 1 + - + 6 

y. 


2y y 2 2y 3 y 


+ e2 ^ " 2 y 2 " 2y 3 + 2y4 + 2 ys) + 


as 


( 2 ) 


(3) 
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The expansion ( 3 ) could also be obtained directly by seeking a perturbation- 
expansion solution of equations (l), that is, by assuming the solution has the 
form 

f = fi(y) + efs(y) + e 2 f 3 (y) + • • • (4) 

and substituting into equations (l), letting e -+ 0 successively to obtain 

the respective problems for £ x , f 3 , etc. From equation (3)j the solution 

apparently is singular as e -> 0 and y -> 0. Also from equation ( 3 ) we see 
that f is analytic for all other y at e - 0, and from equation ( 2 ) we see 

that f is analytic for all other e at y = 0. Thus, the nonuniformity is 

limited to the -neighborhood of the line y = 0, e = 0 in f, y, € space, and 
the function f is, in general, analytic for e = 0 except on that line. 
Finding an expansion valid near y = 0 for small e will be discussed later. 

Let class 2 denote the class of problems for f in which statements la, 

lib, and Ilia apply. An example illustrating this class is example prob lem 2: 


(e + y)f" + (a + l)f’ = 0 , (0 < a < 00 ) 

f(0) = 0 , f(°°) = 1 

4 

The exact solution is 



for which the expansion for small e is 



(5) 

(6) 


(7a) 

(7b) 


The asymptotic solution cannot satisfy the boundary condition at y = 0. Note 
that in this example f is not analytic with respect to e at e = 0 for 
any fixed y unless a is a po stive integer. However, if a = m/n where 
m and n are positive integers, then f is analytic with respect to 
G l = e 1/n at e = 0. Further interesting properties are introduced by the 
variation of the parameter a, to be seen later. 

Subclasses of class 2 have special properties that are of great interest 
in singular perturbation theory. Let class 2a denote the problems of class 2 
for which the singularity in statement lib is an essential singularity at a 
finite location in y space where the boundary condition is lost (without 
loss of generality, say at y = 0) . Consider example problem 2a : 
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( 8 ) 


y^f" + (2y - e)f ’ = 0 
f(0) = 1 , f(oo) = 0 


for which the solution is 


f = 1 


e 


-e/y 


(9) 


This function has an essential singularity at y = 0, and the respective terms 
of its expansion for small e do not satisfy the boundary condition at y = 0, 
although the complete solution does. This solution function can also be said 
to have a boundary layer at y = 00 for large e, which can be seen by using 
the transformation z — l/y, e = i/p , to obtain a problem and solution having 
the same forms for f(z,|3) as the function f(y,e) of example problem pa below. 


Now further denote by class 2b the subclass for which the essential 
singularity is at y = 00 , where a boundary condition cannot be satisfied. 
This subclass is exemplified by "Stokes’ paradox" of viscous flow, and is 
simply illustrated by example problem 2b : 

f" + ef’ = 0 

f(0) = 0 , f(») = 1 



with the solution 


f = 1 


e 


-ey 


= ey 


, € 3 v 3 

21 + 3 '- 


(ID 


The expansion of this solution for small e has a nonuniformity at y = °°; 
the boundary condition there cannot he satisfied by a finite number of terms 
of the expansion for small e. The algebraic approach (in property Ilia) of 
the "inner solution" of the nonuniform region near y = 00 (or near 
z = l/y = 0) toward a finite number of terms of the "outer solution" (the 
expansion of f for small e) as the outer region is approached (l /ey °o) 

is best observed by transforming the problem by z = l/y to the form of 
example problem 2a above. The function (equation (ll)) also has a boundary 
layer for large e at y = 0, as is seen by substituting (3 = l/e to get the 
same form for small (3 as in example prob lem 5a below. 


Denote by class 3 the class of problems for f in which statements Tb, 
lib j and Ilia are true. An illustration of class 3 is provided by example 
problem 2.: 


ef" + (f*) 1 ^ = 0 , (0 < [x < 1) 

f(0) = 0 , f(«>) = 1 


( 12 ) 
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The asymptotic solution found by letting 0 cannot satisfy the boundary 

condition at y = 0 . The exact solution is the same as that for example 
problem 2, given by equation (6) , where: 


a = 


• e = 

[X ’ V n 


e \i-M- 


1-1 J. 


( 13 ) 


Let class 4 be the class of problems for f in which properties lb, lib, 
and Illb apply. Consider example prob lem 4 a : 


e 2 f " + (ky 2 - 6e)( 1 - f) = -2e 
f(0) = f (0) = 0 


( 14 ) 


The exact solution is 


f = 


2y_\ .-y 2 /® 


■ 1/2 


<y/e 


1/2 


; t2 dt 


(15a) 


This function has an essential singularity at e = 0 , but has the asymptotic 
expansion for large y and small e (and the outer expansion) : 


f 



l» 3 e 2 l- 3 ‘ 5 e 3 , 

(Sy 2 ) 2 + (2^) 3 


( 15 b) 


Another interesting problem in this class is example problem 4 b : 

ej^f" + 2y 3 f ' - 2ef = 2y 3 - 4 ey 
f (0) = 0 , f* (0) = 2 

with the solution 



which is regular at y = 0, but has the outer solution, y + e/y. 


(16) 


(IT) 


Let class 5 denote the class of problems for f in which statements lb, 
lie, and Ille apply. A large portion of the quick-transition phenomena of 
mathematical physics, including the classical boundary-layer theory, generally 
fits into this category. A number of examples will be given to illustrate 
various phenomena in this class. 


One of the simplest problems with which to illustrate class 5 is the same 
as equations ( 12 ) but with p = 0. With this example one can observe the 
transition from a problem in class 3 with algebraic approach of the inner 
solution to a finite number of terms of the outer solution (with no essential 
singularity in the solution at e = 0) to one in class 5 with exponential 
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termwise approach of the inner solution to the outer solution and with an 
essential singularity in the solution at e = 0. Consider example problem 5a: 


ef" + f» = 0 
f( 0 ) = 0 , f(°°) = 1 


(18) 


The exact solution is simply 


f = 1 - e 


-y/e 


(19) 


It is instructive to obtain this same result from equation (6 ) , using equa- 
tions (13) and letting p -> 0. In so doing, one needs the limit process: 


lim 
p-*o + y 


1 zUl 

6 \ M. 


= lim ( 1 + ^ 

\i^o\ e 


= lim e 
p-»o 


i (1 * f 


= exp 


- lim 
L |Jr>0 


log(l + py/e) 
M- 


= e 


g 

e 


( 20 ) 


Other simple problems with properties appropriate to class each with 
certain features that will be of interest and useful to study, are the 
following: 

Example problem 5b ‘will be used to illustrate the common occurrence of 
simple first-order displacement: 


€f 1,1 + yf" = 0 

f(o) = f ’ (o) = 0 ; f* (oo) ~ l 


The solution is 


l2ej 


1/2 


f = y erf 

where 

erf | = 2 jc 

The displacement (to be defined) is 


1/2 / 

+ (V^ /2£ - 1 

it 


■1/2 / b -t* 

/ e dt 


= ci/2 ( 2 

it 


5 = e 


1/2 


( 21 ) 


(22a) 


(22b) 


(22c) 


Higher order displacement effects will be illustrated by example 
problem 5c : 
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The solution is 


e ftn + f« = i 

f(0) = f*(0) = 0 ; f’H ~ l + y 


f = y + 4 - - * + 


and the displacement (to be defined) is 


& = e(l-^e + jje 2 - |-e 3 + 


( 23 ) 


( 24 a) 


( 24 b) 


The displacement in a given problem may be complex although both the outer and 
inner solutions are purely real. An example -with a purely imaginary displace- 
ment (to be discussed) is example prob lem 5d : 


ef ,fl + f" = 1 

f (o) = 0 ; f » (o) = 1 ; f ‘ ~ y + o(l) as y ->• «> 


( 25 ) 


with the solution 

f - + e - €e' y/s 

Example problem ^e contains two different displacements (one of which is 
complex) : 

sf't" + f"» = 0 

f(o) = f» (o) = f"(o) = 0 / 

f t ~ y + o(y) as y 00 


(26) 


( 27 ) 


The solution is 

f = ^ - ey + e 2 - e 2 e" y/e (28) 

The cause of more than one displacement and the handling of such problems will 
be discussed. 


A problem in which the asymptotic inner solution is not equivalent to the 
exact solution (to be discussed when the inner solution is defined) is example 
problem 5jT: 


* 

ef" - f = -1 (0 < y < 1) 

f(0) = 2 , f*(l) = 0 


( 29 ) 
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This mathematical model, -which corresponds to a problem in one -dimensional 
heat conduction, was posed by Professor W. C. Reynolds of Stanford University 
(private communication) . A similar problem was discussed by Friedrichs and by 
Van Dyke (see ref. 1, p. 79)- The exact solution to equation (29) is 


f = 1 + . ! + -^ ' y/£ — 

cosh l/e 1/2 1 + e -2/€ 

The asymptotic (inner) solution is 


(y-2) /< 


1 + e‘ 


(30a) 


f ~ 1 + e 


+ exp as e^0in0<y<l 


(30b) 


Internal quick-transition regions and their handling can be illustrated 
by example problem 5g : 

ef" + 2ff = 0 ) 


with the solution 


and by example problem 5h : 


f(-oo) ~ -1 y f(oo) ~ 1 


f = tanh(y/e) 


e(y + e) 2 f" + 2(f - e)[(y + e)f* - f] = 0 
ft (- 00 ) - -1 , f'H ~ 1 


with the solution 


f = (y + e)tanh(y/e) 


and the displacements (to be defined) 


5 =5 = -e 


(3 4a) 


(34b) 


The latter problem can illustrate the calculation of displacement and its role 
in matching concerning a model variable (f*) which is analogous to a variable 
that may overshoot and then relax in a shock wave. 


A class of problems closely related to these previous classes may be 
denoted by class 6 , the problems in which statements lb, lid, and Hid apply. 
Example problem 6: 


ef" + f = 1 


f(0) = f*(0) = 1 


■ 'll llll I II 



with the solution 


f = e 1/a sin (-^ + 1 ( 36 ) 

illustrates this class (see ref. 1, p. 213). Note that f has an essential 
singularity on e= 0 , but not of the exponential type,, and that f f = cos (y/e 1/2 ) 
has no limit as e 0. One does not expect to find physical problems in 
this category. 

A great deal can he learned about the various classes by studying the 
above simple problems. The discussion is limited here because much informa- 
tion is contained in the definition of each class , and the classes and example 
problems -will be discussed later. 

In attempting to solve a given problem by asymptotic or perturbation 
expansion, one may encounter a nonuniformity such as described by one of the 
properties in II. One cannot be certain in advance whether the expansion 
solutions to be determined will be analytic with respect to some form of the 
small parameter. The approach often used in inner and outer expansions is to 
start by assuming a power -series form of the expansion and trying to determine 
the respective terms. In addition, one would like to know which statement in 
III applies, in particular, whether the outer expansion is approached termwise 
algebraically or exponentially by the appropriate inner expansion. If there 
is no discontinuity or lost condition at a boundary (ila) , one can generally 
be sure that there is no quick-transition region and hence no exponential -type 
essential singularity (ilia) . Then, if the terms of the outer expansion are 
singular at y = 0 , the approach will be algebraic (to be explained in the 
subsection The Matching Principle) . If the nonuniformity at e = 0 is due to 
an exponential -type essential singularity, the order of the differential equa- 
tion is reduced, a condition may be lost, and then, if there are no negative 
powers of y in the outer expansion, an appropriate inner expansion can 
approach termwise exponentially to a form of the outer expansion, to be shown. 
Thus, it will be seen (subsection entitled The Matching Principle) that state- 
ment IIIc follows from lie . First- and higher-order boundary -layer problems 
in the fluid-flow theory governed by the Navier-Stokes equations are generally 
of this type. It is in these problems of class 5 where the procedures to be 
outlined below (p. 16 ff) are most advantageously used. 


ESSENTIAL FEATURES OF THE METHOD OF MATCHED ASYMPTOTIC 
EXPANSIONS (INNER AND OUTER EXPANSIONS) 


The purpose of this section is to review briefly the principal ideas of 
the general inner- and outer -expans ion method, which is the primary basis for 
the new approach to be discussed in the next section. For a more complete 
description of the general method the reader should refer to references 1 and 
21 to 25 . 

It is useful to introduce first the most basic principles of asymptotic 
expansions along with the appropriate and more or less conventional notation, 
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for which we may refer to Erd^lyi (ref. 32) . An asymptotic expansion to N 
terms for a function f of a small quantity e and possibly of some indepen- 
dent variables x and y, for example, may have the form 


N 


f(x,y,e) = f n (x,y)cp n (e) + o[cp H (e)] 


n=i 


as € 0 with x,y fixed J 


(37) 


where the "smaller order" symbol o( ) means that, if A = o(B) as e -*■ 0, 

then lim A/B = 0. It should also be specified that cp = 1 and that 
e-»o x 

cp n+1 (e) = o[cp n (e)] as e -*• 0. If f can be expanded (as in eq. (37)) to N 

terms, then one may also write 


M 


f(x,y,e) = ^ f n (x,y)cp n (e) + 0[cp M+1 (e)] I 
n=i I 


as £ 0 with x,y fixed J 


(38) 


where M is any integer from 1 to N - 1 (ref. 32), and the order symbol 0( ) 

means that if A = o(B) as e -* 0 (where B / 0), then lim A/B is bounded. 

6-»0 

An asymptotic expansion for f to N terms may also be written without order 
symbols by using the symbol which represents "asymptotically equal." Thus, 


N 


f(x,y,e) ~ ^ f n (x,y)cp n (e) 


n=i 


as e -> 0 with x,y fixed 


( 39 ) 


The number N may be infinity, but the symbol ~ is still used to indicate 
that f may contain terms of smaller order than any quantity in the sequence 
of cp n as e -*• 0; for example, e _1/€ is smaller than cp n = e n for any 
integer n as e -»■ 0. 


In using the method of matched asymptotic expansions, one usually starts 
with a power series expansion in what appears to be the logical form of the 
small parameter e, and then modifies it as necessary if and when difficulties 
appear. One may also leave the form of the expansions initially unspecified 
so that the expansions can be constructed as the problem proceeds (ref. l) . 

The first step in solving a problem for f(x,y,e) by this method is to 
represent the solution function by the asymptotic expansion (39) where cp 1 = 1 
and N -*■ °°. Equation ( 39 ) is then substituted into the differential equation 
and boundary conditions of the problem. The first-order problem (for fi(x,y)) 
may be extracted by letting e become zero. One then observes that fi(x,y) 
is nonuniform at some values of x and y. Equation ( 39 ) is then denoted as 
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the outer expansion . "Inner variables" are defined by magnifying the physical 
variables in such a way that the inner dependent variables remain uniformly 
valid as e o and as the location of the nonuniformity of the outer problem 
is approached simultaneously. The inner variables are thus required to be 
0(l) as e -> o. If the nonuniformity is at y = o when e = o, then the 
independent variable y is magnified by some function of € as 

y = cTi(e)Y (40) 


The dependent variable f should be similarly magnified: 

f(x,y,e) = o- 2 (e)F(x,Y,e) (41) 

•where the inner dependent variable F is o(l) as € o with x and Y fixed. 
The functions a i and cr 2 must be determined according to the requirements of 
the specific problem upon substituting equations (40) and (41) into the 
governing equations and boundary conditions. 


An "inner" asymptotic expansion for F, valid in the inner region, is 
also left unspecified a priori and may be written as 


F(x,Y,e) ~ ^ F n (x,Y)<i> n (e) 

n =i i 

as e 0 with x,Y fixed 


(42) 


where $>! = 1 and 0 n+1 = o(o n ) as s -»• 0. As discussed above, one may try 
letting the sequence of 0 n (e) be a power sequence of what appears to be the 
appropriate form of e, until difficulty is encountered, but, in general, the 
form of each $ n (e) ™.s.t be determined as one proceeds with the problem. 


The solutions may be obtained for the respective terms of the inner and 
outer expansions when sufficient boundary conditions for each of the problems 
are known. Sufficient boundary conditions are provided by a procedure for 
matching the inner and outer expansions to each other according to the princi- 
ples set forth by Kaplun (ref. 21) (see discussion by Lagerstrom, ref. 25). 

In essence, the outer expansion of the inner expansion matches the inner 
expansion of the outer expansion to appropriate orders for the respective 
terms. The matching principle also can often be used to determine the inner 
and outer asymptotic sequences of functions of e (ref. l) . 


The matching principle is conveniently stated with the following notation 
Let f mo be the m term outer expansion for f, that is, 

f mo = fi + cp2( e ) f 2 + • . . + cp m (e)f m (43) 

and let f ni be the product of cr 2 (e) times the n term expansion of F, 
that is, 

f ni = CT 2 [ Fi + 0 2 (e)F 2 + . . . + $ n (e)F n ] (44) 
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If we now substitute inner variables into the outer expansion (43),. we have 
"the m term outer expansion in inner variables / 1 denoted as (fmo)i. if we 
substitute outer variables into the inner expansion (44), we have "the n 
term inner expansion in outer variables," denoted as (f ni )°. Then the first 
n terms of (f 1110 ) 1 constitute the /'n term inner expansion of the ^m term 
outer expansion," denoted as (f^ 0 ) 311 , and the first m terms of (f ni )° con- 
stitute the "m term outer expansion of the n term inner expansion," denoted 
as (f ni ) mo . The principle that determines the matching between the inner and 
outer expansions is then 


( f mo)ni _ ( f ni)mo ( 45 ) 

The use of intermediate expansions and intermediate matching is also often 
necessary or convenient (see ref. 1 or 23). 

After the inner and outer problems have been solved it may be desired to 
combine the inner and outer solutions into one solution called the "composite 
solution," which approaches the outer solution in the outer limit and 
approaches the inner solution in the inner limit. The composite solution is 
therefore uniformly valid over the entire region of interest. In the above 
notation, the "additive rule" for forming the composite solution is 

fC _ f ni + f mo _ ( f mo)ni (l^g) 

METHOD OF ASYMPTOTIC EXPANSIONS WITH DISPLACEMENT VARIABLES 


In this section will be outlined a new approach to obtaining inner and 
outer asymptotic -expansion solutions in certain classes of problems that have 
been formulated in terms of differential equations and boundary conditions and 
that contain a small parameter. The essential features distinguishing this 
approach from the conventional inner- and outer -expansion method are: 

(a) a specific definition of outer and inner functions as formal sums of 
series that converge in some region for sufficiently small values of the 
expansion parameter and which are asymptotic to the exact solution; 

(b) the use of displacement variables both in defining the outer solution 
function and in matching the inner solution to the outer solution; 

(c) the development and statement of a matching principle that explicitly 
retains all displacement effects in matching inner solutions to outer solu- 
tions in a class of problems; and 

(d) a statement of certain rules often found to be useful in obtaining 
the appropriate forms of the expansions and magnifying factors. Application 
of the approach formulated here is restricted to certain classes of problems 
(to be discussed) . 
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To show how this method takes advantage of the special properties of 
class 5 (as defined in the section "Characteristics of Several Classes of 
Singular Perturbation Problems") and to show the distinctions from other 
classes, the solution of problems in classes 1 and 5 will be discussed. 


Analytic Inner and Outer Functions as 
Sums of Assumed Power Series 


Sums of asymptotic series and of analytic expansions .- Consider a 
function f (x,y,e) with the asymptotic expansion : 




00 

c P n ( € ) f n( x >y) 

n=i 


as e 0 


( 47 ) 


where 


cpi = 1 and cp n (e) = o[cp n . 1 (e)] as e 0 
and where, for each cp n (e), there is some m > 0 such that 

e m = o[cp n (e) ] as e - 0 


(W) 


The series on the right side of (47) is an asymptotic series . The series 
may be convergent or divergent (ref. 32, p. 12). Denote by D*(x,y) the 
domain of all values of x,y for which the series converges for sufficiently 
small e > o. Define f*(x,y) (in D*) as the sum of the infinite series of 
functions cp n (e)f n (x,y) , excluding all parts of f that are o[cp n (e)] as 
e o for all n (e.g., terms that are exponentially small if the cp n are 
powers of e, cp n = e n ) : 

L1 oo 


f*(x,y,e) 


c P n ( e ) f n( x >y) in D *( x >y) 


( 49 ) 


n=i 


Although the expansion (47) does not uniquely determine its asymptotic sum f 
(ref. 32, p. 14) because f may contain terms of smaller order as e o 
than any cp n (e) of the asymptotic sequence, the sum of the series, f*, is here 
defined to contain only the terms cp n f n an -d so ^- s uniquely defined by 
equation (49) j if the series converges. 


From this point on, only continuous functions that possess asymptotic 
expansions in series of powers of some function e c of the small parameter 
will be considered (except in the subsection "Extension of the Method to a 
Class of Problems With Outer and Inner Functions Nonanalytic at e = 0," 
pp. 33-36); thus it is assumed that f has the asymptotic expansion 


f(x,y,e) 


OO 

^ £c _1 f n (x,y) 


n=i 


as e -* 0 


(50a) 
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where 


e c = e c( € ) 0 as e -*• 0 (50b) 

and where e c (e) is determined in some way. (The subscript c does not yet 
have any special meaning here, but is convenient for later use.) For all 
x,y in D*(x,y), defined now as the domain where the series in (50a) con- 
verges for sufficiently small e > 0, the sum of the series excluding terms of 
f that are of order smaller than e c m for any m is an analytic function of 
e c at e = 0, denoted as 

00 

f*(x,y,e) = Y €g -1 f n (x,y) in D*(x,y) (50c) 

n=i 

(from the theory of analytic functions of a complex variable) . Conversely, 
if, in some domain D*(x,y) , f* is analytic with respect to e c at € = 0, 
the series (50c) converges for sufficiently small e there and is identical 
to the Maclaurin series in € c . 

Denote the difference between f and f* by q*(x,y,e): 

f(x,y,e) = f*(x,y,e) + q*(x,y,e) (50d) 

Since q.* contains only the parts of f that are smaller than any power of 
e c as e -*■ 0 (by the definition of f *) , then 

q.*(x,y, e) ~ 0 + exp as e -* 0 (50e) 

where "exp as e -»■ 0" indicates terms that are exponentially small as e -» 0 
with the other arguments (x,y) fixed. 

The analytic function f*(x,y,e) is defined in a nonvanishing domain 
D*(x,y) if, as assumed, f has the asymptotic power-series form, (50a), and 
if also at some boundary, say y^ (x^) , a boundary condition is applied to f 
that is an analytic function of e c at e = 0: 

CO 

f ( x b>yb> 6 ) = f *( x b>yb,e) = Y eJ^fnCxjj^b) (50f) 

n=i 

(The above statement is easily proved with the comparison test for convergence, 
with the geometric series as majorant.) 

An example that conveniently illustrates the above discussion is the 
function 
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_ / x r y dt -y/e 

f(y ’ e) * J mt - 8 + 


-X/€ 

e > 


(y > 0) (51a) 


= f*(y,e) + q.*(y>e) > (0 < y < °°) 


( 51 b) 


where 


f*(y,e) - J V f4f - £ (-l) n nU n [l-(j- + y + |r + ---+Sr) e 


1 + et 

° n=o 


z~\ a ~y 


and 


l*(y,e) = e' y/G + e" l ' e 


( 51 c) 

( 51 d) 


The function f(y,e) is not analytic with respect to € at e = 0 , hut its 
asymptotic expansion for small e is the asymptotic power series in (51c). 

The series ( 51 c) converges for |e| < | l/y] , and so f* is analytic with 
respect to e at e = 0 for all finite y, and D*(y) is then 0 < y < 00 • 

For y = 00 , the function (l + et) -1 e~^ dt becomes nonanalytic at e = 0 , 

since the radius of convergence about e = 0 in the complex e plane vanishes, 
and the asymptotic expansion 


f(“,e) 



e' -1 dt 
1 + et 


1 


lie + 21 e 2 


3 *e 3 + 


( 51 e) 


diverges (cf. ref. 32, pp. 1-3) . It is pertinent to note here that, in this 
typical example, f* is not a good approximation to f very near y = 0 for 
small e, since at y = 0 , f = 1 + e“ 1/€ and f* = 0. Now let y = eY and 
consider the same function f(y,e) as in equation (51a) but with y replaced 
by eY : 


f(eY,e) = 


‘eY e-t dt 


1 + et 


-Y 

+ e + e 


■i/e 


(Y > 0) 


(52a) 


This function of Y and e, which is identically equal to f(y,e), is also 
nonanalytic with respect to e at e = 0, but has the asymptotic expansion 


f(eY,e) ~ e“ Y + ^ Aj(Y) e 2j+1 + ^ 


B i (Y)e 2j+2 


(52b) 


J=o 


0=o 


where : 
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(2j + 1 - kJT 


9 


(52c) 



k=o 


Bj(y) 



(-j) k k , .Y gJ+g " k 

(2j + 2 - k)» 


The two series in (52b) converge for all |e 2 | < |l/yj. Thus, f(eY,e) can he 
written, by analogy with equation (50d), as 


f ( eY , e ) = P(Y,e) + Q(Y,e) , (0 < Y < °°) (53&) 

where F(Y,e) is the analytic function in e at e = 0 which is the sum of 
the power series in e, containing all parts of f not exponentially small as 
e -*■ 0 with Y fixed: 

F(Y,e) = J (53*>) 

o 

and where 

Q(Y,e) = e 1/6 ~ 0 + exp as e -» 0 with Y fixed ( 53 c ) 


The domain of definition of F(Y,e) (like f*(y,e)) is (0 < Y < 00 ) . The 
function F(Y,e) is a good approximation to f(y,e) near y = 0 for small e. 

In another example (cf. example problem 5gj eqs. (31)), 

f = tanh(x/e) (54a) 


has the asymptotic expansions: 


f ~ 1 + exp as 

e -*• 0 

(x > 0) 


f ~ -1 + exp as 

e -*• 0 

(x < 0) j 


In this case, two asymptotic analytic functions 

f* are 

needed: 

f*+ = 1 , x > 0 ; 

f*- = 

-1 , x < 0 


(54b) 


(54c) 


since the asymptotic solution is discontinuous at x = 0. This is an example 
of Stokes’ phenomenon (ref. 4). For this case, the analytic function of e 
that is uniformly valid in the region of the nonuniformity of f(x,e) is 

F(x) = tanh X (54d) 

where X = x/e is held fixed as e -*• 0. 

In solving for an unknown function f in a singular perturbation 
problem, it is not known for certain whether asymptotic expansions of the type 
(50a) can properly represent the function. If one assumes that such an expan- 
sion does represent the function, then one can speak of the analytic sum f* 


20 



of the assumed power series (or, more "briefly, assumed - analytic function ) . 

If the form (50a) is valid as assumed, then the function f* defined "by 
equation ( 50 c) exists at least in some neighborhood D* of a boundary where a 
condition is applied to f that is analytic with respect to e c at e = 0. 

In the approach to finding inner and outer expansions used here, assumed - 
analytic functions (like f* and F in equations ( 51 ) to (53)) 'will be used for 
both the outer and magnified inner functions. The outer function is assumed 
to be a power series in some parameter 6 q = € o( £ ) the magnified inner 

function to be a power series in = e^(e). These power series are assumed 
to represent the solution asymptotically as e 0 and are analytic functions 
of e Q and e^, respectively, at e = 0 in the respective regions where the 
series converge. If consistent solutions cannot be determined for the respec- 
tive terms of an expansion, then the assumed-power-series function is not 
applicable, that is, cannot represent the exact solution function asymptoti- 
cally. If the various terms in the assumed expansion can be determined and 
are consistent and noncontradictory, then it is assumed that the function f* 
given by the expansion represents the function f asymptotically. This is 
equivalent to an assumption of uniqueness of the asymptotic expansion, to. the 
order in € c considered . The asymptotic expansion of a function is known to 
be unique Tor a given asymptotic sequence (ref. p. 13 )* If* in a given 

problem, there is reason to doubt that outer and magnified inner functions 
that are analytic in some e c (e) at e = 0 can be defined even though the 
first few terms of a power series in e c are found, then a slightly different 
point of view may be adopted (see subsection "Extension of the Method . . 
pp. 33 - 36 ), but the terms already found remain valid to the order considered. 

In many problems either the outer function or the inner function or both 
may be equivalent to the exact solution function. For instance, example 
problem 1 will have both outer and inner functions equivalent to the exact 
solution function; in class 5j however, the outer solution cannot be equiva- 
lent to the exact solution, but in all the example problems of class 5 except 
5f the inner solution function is the exact solution. In example problem 5f 
neither the outer function nor the inner is the exact solution. 

Definitions of inner and outer functions . - The definitions of the outer 
and inner functions to be used in the present approach are contained in the 
following paragraphs. 

Assume that f(x,y,e) can be separated into two parts, one part (f°) 
which is analytic in some region with respect to some form e 0 (e) at e = 0 
and the other (q°) decreasing to zero faster than any power of e Q : 

f(x,y,e) = f°(x,y,e) + q°(x,y,e) (55a) 


where 

f° is analytic with respect to e Q = e Q (e) at e = 0 in some D°(x,y) 

q° ~ 0 + exp as e ^ 0 with x,y fixed J 

(55b) 
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Note that only terms contributing to f° will appear in any asymptotic power- 
series expansion of f as e 0 with x,y fixed. Except near some nonuni- 
formity of f at some (x,y), f° is then a valid representation of f and is 
defined as the outer function (if, in fact,, f has an asymptotic power series). 

If the order of the differential equation is not reduced and all the 
boundary conditions can be satisfied, one can treat the problem as class 1. 

In this class there are no exponentially small terms in the exact solution, 
so q° = 0 and f° and its expansion for small c are equivalent to f : 

f ° = f , (class 1) (56) 

The outer expansion is then determined as the power series in e in the form 
of, say, equation (50c) which, by direct substitution, is made to satisfy the 
differential equation and all the boundary conditions on f(x,y,e). This 
expansion may then be not uniformly valid near some x,y because individual 
terms are singular there . 

In boundary -layer -type problems (class 5) the order of the differential 
equation is reduced upon letting €->()• Then the limit solution obtained is 
found to be discontinuous; that is, either it does not satisfy a condition of 
continuity of f (or one of its derivatives) at an internal location (called a 
Stokes line; cf. ref. 4), or it cannot satisfy (continuously) a boundary con- 
dition required of f or one of its derivatives. The discontinuous behavior 
of the limit function indicates the presence of terms in f that vary rapidly 
(exponentially) near the location where the discontinuity is encountered. The 
discontinuity itself is an asymptotic representation of a rapid variation of 
the dependent variable over a very small range of the independent variable 
(say over y near y = 0) . 

An outer function f° that corresponds to the nonuniformity at y = 0 
in a problem of class 5 may be defined by equations (55); where q° contains 
the exponentially decaying quantities (as e 0) including those that cause 
the rapid variation of f and its derivatives near y = 0. (in the case of 
an internal discontinuity, the function f° in equations (55) is given by two 
different analytic functions, f° + and f°“, one valid on each side of the dis- 
continuity.) Since q° and its derivatives are exponentially small everywhere 
except very near y = 0 when e is small, f° is a valid asymptotic function 
representing f everywhere except very near y = 0. Then f° must satisfy, 
to within exponentially small amounts , the differential equation and all 
boundary conditions of f and its derivatives except those at y = 0 . Further 
conditions (at or near y = 0) that determine f° are discussed below. 

Denote by w(x,y,e) the particular dependent variable, or derivative of 
f, whose limit function in the outer problem is discontinuous; also define 
w° as the corresponding function in the problem for f° and let 

w(x,y,e) = w°(x,y,e) + q. w (x,y,e) (57) 

where ~ 0 + exp as e -* 0. (in example problem 5b the outer limit func- 
tion is fi(y) = y; the derivative that cannot satisfy a condition at y = 0 
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is n(y) = i, so , in that problem, w = f f (y,e) 9 ) At y = 0, neither the 
complete function f° nor any of its derivatives can generally satisfy 
exactly the same conditions as are satisfied by f and its derivatives 
because of the significant terms in q° and its derivatives at y = 0. Con- 
sider example problem 5b, for which the solution f is given by equation (22a) 
and f f = erf(y/(/2e), and in which f and its derivatives are represented 
asymptotically to within exponentially small quantities by the outer functions 

f°(y) = y - e 1/2 v/2/^ , f°* = 1 , f°" = 0 ( 58 a) 

At y = 0 ; 

f = 0 , f« = 0 , f" = s“ 1/2 n/2 , f m = 0, . . . 

f° = -e 1/2 N/2^T , f°* = 1 , f°" = 0 , f° ,n »0,... 4 

( 58 b) 

Note that the deviation of w° from w = f 1 at y = 0 is o(l) as e 0, 
whereas the deviations of the derivatives of lower order than w° are o(l) 
as e -> 0. One possible way to regard this is that, to first order, only the 
condition on w = f f must be omitted from the outer problem, but that, to 
order e, the condition on f is not satisfied. But then there would not be 
sufficient conditions to determine f° to second order. A more useful 
approach is to omit only the condition on w from the outer problem and to 
require f° to satisfy the same condition as f (i.e., f = 0 in this problem) 
in some way, to ma ke f° determined. But we see that f° vanishes at 
y = 5 = € i/s J 2 /it rather than at y = 0. Thus, we require the conditions on 
all lower order derivatives of f than w to be_ satisfied in the prob lem for 
f° , but allow those conditions each to be satisfied at an appropriate displace- 
ment surface , y = 5 . In general, there is a displacement associated with each 
condition applied to f° near y = 0. Corresponding to each condition 

y = 0 , d m f(x,y,e)/dy m = (59a) 

for all lower order derivatives than w, tiie displacement & for the function 
S IT1 f 0 (x,y,e) /3y m is therefore the root 6 of 

d m f°(x,S,e)/dy m = % (59b) 

Since 5 is small, the conditions can usually be transferred by Taylor T s 
series to y = 0. The displacements are unknown a priori in solving a problem, 
but are determined as part of the solution (to be shown) . One must be careful 
to note that even in problems involving only real variables with purely real 
solutions, the displacements may be complex. In example problem 5e, for which 
w = f"(yj€) and for which two conditions must be applied in the problem for 
f° near y = 0, the condition f° = 0 is satisfied at y = 5 a = e(l ± i) and 
f°* = 0 is satisfied at y = ^ = e. (Complex displacements in general need 
not be considered unless a difficulty is encountered not otherwise explainable 
in matching inner solutions to outer solutions or in determining some unknown 
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quantity.) With the above requirement that all conditions on dependent 
functions of lower order than w he applied at displacement surfaces , the 
outer function f° is determined in terms of unknown displacements in class 5- 

Although the individual terms of the expansion of f° will he determined 
from reduced-order equations, the function f° is still determined hy the 
complete differential equation. The condition on w at y = 0 is replaced hy 
the condition that there he no boundary layer, or rapid transition, in the 
function f° near y = 0. This insures a unique solution for f°. 

In the inner region of the nonuniformity (the vicinity of the location 
where the terms of the outer expansion become infinite in class 1 or the 
vicinity of the discontinuity of the outer solution in class 5)> an asymptotic 
function representing the solution function is to he found hy magnifying both 
the independent and dependent variables in such a way that the solution 
remains uniformly valid as e -*■ 0. Let the magnified independent variable he 

Y = y/oi(e) (60a) 

and let F = f 1 /cr2(e) denote the assumed -analytic function (cf. discussion of 
f* above) representing (at least asymptotically) the magnified solution 
function, f/o2(e)> in the inner region, that is, 

f (x,y,e) = f 1 (x,Y, e) + q 1 (x,Y,e) (60b) 


where 


F(x,Y,e) = f i (x,Y,e)/a2(e) is assumed to he 
analytic with respect to some quantity 6j_ = €j_(s) 
at e = 0 in some D^(x,Y) ; and 
q 1 ~ 0 + exp as e -*■ 0 with x,Y fixed. 


(60c) 


All inner dependent variables (F,BF/dY, etc.) are considered to he 0(l) as 
e -» 0 with Y held fixed. Substitution of equation (60b) into the differen- 
tial equation of the problem for f and into the boundary conditions on f 
and its derivatives at y = 0 requires that fi(x,Y,e) satisfy those equations 
and conditions to within exponentially small amounts as e -> 0 with x,Y 
fixed. The problem for F(x,Y,e) is then determined if appropriate forms for 
ax(e) and p 2 (e) are found and if sufficient outer boundary conditions are pro- 
vided by appropriate matching of the inner function to the outer function. 

Some particularly useful matching relationships can be derived from the defini- 
tions of the outer and inner functions (to be shown below) . 


The functions e D (e), ei(e), Oi(e), and 02 (e) must Le determined 
according to the specific problem. Some rules that are often useful for 
deter mi ning 01 and 02 will be discussed later. When oi and 02 are known, 
and the inner problem is set up in terms of x,Y,F, and e, one can assume e Q 
and ej_ to be the forms of the small parameter occurring explicitly in the 
completely - specified outer and inner problems (including the conditions to be 


24 



applied at the displacement surfaces in the outer problem in class 5 and the 
conditions to be supplied by the matching principle in the inner problem) . 

For the most convenient use of the matching principle, both the outer and 
inner problems may be solved in terms of a single small parameter e c = e c (e) 
with respect to which both e 0 and are analytic. For example, if 

€ Q = e 1/2 and = e 1/3 , then the common form of the small parameter to use as 
the expansion parameter is e c = e 1/s . The assumed power series are then 

f°(x,y,e) = fi(x,y) + e c f 2 (x^y) + e c 2 f 3 (x,y) + . . . (6l) 

F(x,Y,e) = Fi(x,Y) + e c F 2 (x,Y) + e c 2 F 3 (x,Y) + . . . ( 62 ) 

which will be substituted into the inner and outer problems so that the 
respective terms can be evaluated. In the subsection "Extension of the 
Method . . .," (pp. 33-36) the requirement of analyticity of f° and F at 
e c = 0 is relaxed, but the terms of equations (6l) and ( 62 ) remain valid if 
determined. 


Magnifying Factors 

The magnifying factors a i and a 2 in equations (60) define the relation- 
ship between the inner functions F(x,Y,e) and fi. The order of magnitude of 
cfi represents the order of the width (measured in the units of y) of the 
region of nonuniformity, and a 2 represents the order of magnitude of f in 
that region as e -> 0. 

For class 1, a solution that is valid near the singularity at y = 0 as 
e -» 0 must be found by defining an inner function that is, in general, ana- 
lytic as y -► 0 and e 0 according to equations (60) . In the problems of 
class 1, the outer expansion (expansion of f°) is the same as the formal 
perturbation expansion of f (x,y,e) . The proper magnification factors a ± 
and 02 are therefore determined by comparing successive terms in the outer 
expansion, observing how the nonuniformity "grows" (for every order of e 
added to the expansion, a certain higher -order infinity is encountered; see 
ref. l) , then observing the order of the dependent variable as the singularity 
is approached, and defining a i and a 2 so that Y and F are of order unity 
in the inner region (as e and y 0) . The terms of the inner expansion con- 
tain constants of integration that must be evaluated by matching the inner and 
outer expansions. The matching principle is introduced in a later section, 
but it may be noted that for problems of class 1, the matching process is 
equivalent to that used in the conventional inner and outer expansion method 
except for the particular use of the common parameter e c here. 

In problems of class 5* some simple rules have been found to be useful in 
determining the magnifying factors a ± and 0 2 . If these rules lead to a use- 
ful result without contradictions in a given problem, then the choices made 
for 0i and c 2 are valid. The first rule is: If a function denoted by w 

is discontinuous (but bounded) in the outer problem (e.g., if w° cannot 
satisfy a boundary condition required of w as e -> 0 with y held fixed) , 
then require 
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wi = 0(l) as e -*• 0 with Y = y /a, held fixed (63) 

The function w may he f itself or some higher order y derivative of f . 
For example (as in example problem 5b), suppose w° = df°/dy cannot satisfy 
the boundary condition required of df/dy at y = 0. Then require that 

wi = df^/dy = o(l) as e -*• 0 with Y fixed 

Since all the magnified inner functions (F and its Y derivatives in 
example 5b) are required to be 0(l) as e 0 , we can then write 


dF _ oh df£ 

dY 02 dy 


0 ( 1 ) 


as € -* 0 


For example problem 5b , the result is that 


tfi/02 - 0(1) 

and we may let a(e) = a i(e) = 02(e). 
is as follows (Refer to fig. 1.): The 

y 



as e -* 0 

The reasoning (for a general problem) 
discontinuity of w° is an asymptotic 
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Figure 1.- Illustration of variable w that has internal discontinuity in the outer problem. 

representation of a rapid variation of w. To study the variation of w in 
the inner region represented by that discontinuity, one must magnify the y 
scale (Y = y/ai) but keep the w scale fixed as € 0. This is so because, 

if the order of magnitude in the w scale were changed as e 0, no varia- 
tion of w could be observed (on the right side of fig. l) because all the 
values of w^ being considered would either become infinite or vanish as 
€ -> 0. Hence, in the inner solution we need w^ = o(l) as e-»0. The type 
of function for which w is discontinuous at a boundary can also be referred 
to for the above discussion. In the same scheme as figure 1, that function 
would be as shown in figure 2. There remains one unknown magnifying factor, 
which can be found by the second rule, the princ iple of least degeneracy 
(attributable to Van Dyke (ref. 1) in the general method of matched asymptotic 
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Figure 2.~ Xlius'bi'a.'fcion of verieijls w ■fcha.’fc has lost hound-axy cond.i"fcion in "the outer probiem. 

expansions). Equations ( 60 ) are substituted into the differential equation 
and the result is examined for the limiting process e->0. As € -* 0, the 
equation must not degenerate into a lower order differential equation and the 
most significant terms must not be lost in the first-order problem. (For 
example problem 5 b , then, a = e 1 / 2 .) 

In problems of class 5, the outer problem is often influenced by the 
inner problem, that is, there are displacement effects. Hence, e Q and 
cannot be determined before the inner and outer problems are completely speci- 
fied, including the inner boundary conditions that are not lost in the outer 
problem, to be applied at the displacement surfaces, and including the outer 
boundary conditions to be supplied to the inner problem by the matching prin- 
ciple. The independent displacement variable y (to be defined) will be seen 
to contain e, so that e then occurs in the boundary conditions of the outer 
problem and will also influence the inner problem through the matching process. 
The choice of e c therefore depends on the form of the displacement 
variable y. 


Displacement Variables 

The influence of the nonuniform region on the outer solution is known as 
the displacement effect. The most common example occurs in fluid boundary 
layers. The flow around a body in a uniform stream of fluid would be an 
inviscid flow if the body were frictionless and if there were no other 
source of viscous effects. If, however, the fluid flowing at high Reynolds 
number (Re -> °°) adheres to the surface, most of the flow is inviscid but the 
viscous effects very near the surface produce the quick-transition region or 
boundary layer. The outer inviscid flow is then displaced a small distance 
from the body in relation to the corresponding flow over a frictionless body. 
There is a surface, which we denote as the displacement surface , that would be 
the equivalent frictionless body for the outer displaced flow. (Refer to 
Lighthill, ref. 33 , for a discussion of displacement; also see Mangier, 
ref. 34-.) 
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Any mathematical problem that contains a displacement effect (cf . example 
problems 5b and 5c) can be thought of as having a displacement surface (or 
more than one, depending on the number of inner boundary conditions). Let y 
be the coordinate normal to and measured from the actual surface where the 
nonuniformity (discontinuity or lost boundary condition) occurs and y (the 
displacement variable) be the coordinate measured from the displacement sur- 
face, whose location is unknown. The displacement is then y - y and is, in 
general, a function of the other independent variables (say x and e) : 

S = 0g(e)A(x,e) = y - y { 6 k) 

where cTg(e) is the order of magnitude of 5 as e -*■ 0. If the displacement 
is real, it is generally of the same order as the region of nonuniformity 
(o(ai) as e -* 0), that is 

cXg = 0i for 5 real (65) 

(However, in the special case where 5 is complex, it may have a different 
order of magnitude than 01; in example problem 5b, 01 = e but 5 = o(ie 1/£ ).) 

Since the outer function f° contains no exponentially small terms, the 
function b(x,e), as the root of equation (59b), also will contain no exponen- 
tially small terms. It is then compatible with the assumption that f° is 
analytic with respect to some e c to assume A is also analytic with respect 
to e c . Thus, 

A(x,e) = Ai(x) + e c A 2 (x) + e c 2 A 3 (x) + . . . (66) 

(where, indeed, this requirement is to be accounted for in determining e c ). 
The terms of A are to be determined asymptotically by the inner solution 
F(x,Y,e). 

Useful relationships between the outer displacement variables and the 
inner variables are: 

y = 01Y - 0g[Ai(x) + e c A 2 (x) + e c 2 A 3 (x) + . . .] (67a) 

x = x (67b) 

Y = y/01 + (a & /0i) [A x (x) + e c A 2 (x) + e c 2 A 3 (x) + . . .] (67c) 

For use of displacement variables in matching, define an outer displace- 
ment function f°, which is a function of x, y, and e and is given by 

f°(x,y,e) = f°(x,y,e) (68) 

and its expansion by 

f° = f i(x,y) + e c f 2 (x,y) + e c 2 f 3 (x,y) + . . . (69) 

Note that equations (68) and { 6 k) together, along with x = x, simply 
represent a transformation of variables, that the power -series expansion (69) 
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is equivalent to ( 6 l) (because of eq. ( 68 )), “but that generally the f n are 
not the same functions of x,y as the f n are functions of x,y (except, of 
course, for the first term, f x ) . Naturally, n terms of the expansion of f° 
agree with n terms of f° to the appropriate order (i.e., to C^eg" 1 )). 

Note also that a set of displacement variables corresponding to each inner 
boundary condition must be applied to f°. In a large class of problems, how- 
ever, there is only one displacement. In that^case, the entire outer problem 
may be transformed to displacement variables (f°,x, and y) before solving it, 
or it_ may be solved in terms of x and y, then transformed to displacement 
variables (see application to blunt -body fluid flow problem below). 


The Matching Principle 

Outer boundary conditions to be satisfied by the respective terms of the 
inner expansion mast be obtained by matching the expansion of the inner func- 
tion term by term to the expansion of the outer function. The matching 
referred to here is an analytical matching to obtain rigorously the form of 
each term in an expansion, in contrast to frequently used approximate methods 
whereby functions or expansions are "matched" (or "patched") at some arbitrary 
point (cf. discussion by Van Dyke, ref. 1 , p. 89^ and Chang, ref. 35 ; 
pp. 819-820) . The matching principle (primarily due to Kaplun) in the conven- 
tional inner and outer expansion method is expressed by equation ( 45 ). The 
displacement variables to be used here are not employed there and the displace- 
ment effects do not enter into the process of matching. The following para- 
graphs discuss several proposed forms of the matching principle to apply in 
the present approach to constructing inner and outer expansions. These pro- 
posed forms are similar to Kaplun 1 s proposition. It has been found that, in 
problems of class 5,? a form of the matching principle can be applied that 
retains explicitly all displacement effects. The nature of the displacement 
effects, as well as their role in the relationship between the inner and 
outer solutions, is then more easily and completely determined than in the 
conventional method. 

For this discussion of matching, let superscript i denote a function 
that has been transformed from outer to inner (magnified) variables and super- 
script o denote a function transformed from inner to outer variables, in 
addition to the previous use of f° and fi to denote outer and inner func- 
tions. Then (f 0 ) 1 is the entire outer function transformed to inner vari- 
ables. A combination of equation ( 55 a) written in inner variables and 
equation (60b) gives 

f 1 = (f 0 ) 1 + (q 0 ) 1 - q 1 ( 70 ) 

The difference (q 0 ) 1 - q 1 consists only of the terms of q° that, when 
transformed to inner variables, become of larger order than o(e rri ), where m 
is some finite number, as e 0; but those terms are then exponentially small 

as Y -» °° . (in the example of equations (51) to ( 53 ) j q 1 - (q 0 ) 1 = -e‘^.) 

Thus, 

f i „ ( f O)i + exp as Y -* 00 (71) 
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A direct corollary of (71) is that, to order o^e S'" 1 in inner variables. 


f ni _ (fO)ni + exp 


c 

as Y 


( 72 ) 


where superscript ni indicates n terms of an expansion written in inner 
variables. _ Note that for class 1 qi = q° = 0 and there is no difference 
between f 1 and (f 0 )^. (in class 1 the only difference is in the expansions.) 
The asymptotic relations (71) and (72) are valid in regard to any function f 
for which f° and fi can be defined by equations (55) and (60). Equa- 
tions (71) and (72) are useful in themselves only if the complete unexpanded 
outer function (or its form) is known. However, generally useful forms of the 
asymptotic relation for the termwise matching of expansion solutions can be 
deduced directly from the form (72) (as seen below) . 


In the present matching procedure, when an outer expansion is to be 
matched to an inner expansion, the outer will always be written in terms of 
the displacement variables x,y. The superscript "no” will indicate n terms 
of the expansion of f°(x,y,e) (see eqs. (68) and (69)). Thus, the n term 
expansions to be matched are: 

f no = ?i(x,y) + e c f 2 (x,y) + e c 2 f 3 (x,y) + . . . + e£ _1 f n (x,y) (73) 

f nl = 0 2 (e)[F!(x,Y) + e c F 2 (x,Y) + . . . + eg _1 F n (x,Y) ] ( 74 ) 

If the transformation by equations (67a) and (67b) is used in ( 73 ) j and if the 
transformation by equations (67b) and (67c) is substituted into (74), one 
obtains (f 110 ) 1 and (f ni )°, respectively. If these expansions are then trun- 
cated to n terms (i.e., to order eS** 1 in outer variables and to order 
o'2e£“ 1 in inner variables), one has (f no ) ni and (f ni ) no . 

If the first n terms (to order e^" 1 ) of an expansion of f° contain 
all the terms of f° that would appear in the first n terms in inner vari- 
ables, (f°) ni , then the quantity f° in (72) can be replaced by f no . This 
condition is apparently satisfied in problems of class 5. it all of the 
following (to be discussed) are true: 


(a) the outer expansion is first written in ^ 

displacement variables (eq. (73) ) 5 

(b) there is only one displacement in the function 
being matched; and 

(c) 5 = o(cji) as € -► 0 . 


(75) 


Then 

f n ^ ~ (f 110 ) 11 ^ + exp as Y -* c 0 , (class 5) (76) 

In class 5 (as defined in pp. 4 - 13 ) , there are no negative powers of y in 
the expansions of the f for small y. Then the terms of f° that are 
0(e^ -1 )^ when transformed to inner variables, will become at least as small as 
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order for each n. Then, if ( 75 a,b,c) are satisfied, the terms of 

the displacement in the outer solution that will affect the nth-order inner 
problem are not lost upon truncation of the outer expansion at n terms and 
the terms of the displacement that affect the nth inner problem are also 
carried along in the transformation of f no by equation (67a) to inner vari- 
ables. (These statements are easily checked in any given problem.) The 
respective terms of the expansion of the displacement can then be evaluated in 
the inner solution to satisfy (76) exponentially. In example problem , 
f 0 (y,e) = y - e 1/2 Ai and f 0 = y = f xo . Thus, f x ° contains the first-order 
displacement (the entire displacement in this case). Transformation to inner 
variables then gives (f xo ) xx = e 1/2 (Y - Ax), which carries along the first- 
order displacement, to be evaluated by compatibility with the first inner 
solution. If ( 75 b) is not satisfied, that is, if f° contains more than one 
displacement in a given problem (because of more than one displaced boundary 
condition to be satisfied by the outer function) , then all the nth order 
displacements cannot in general be carried through the above process without 
some loss. However, if the terms of only one displacement (expanded in g c ) 
appear in an expansion of some derivative or function, then the matching rule 
can be used for that function using the appropriate displacement variables. 
(This occurs in example problem 5 e.) If ( 75 c) is not satisfied (i.e., if 6 
is larger than 0(ai) as e ^ 0) , which is possible if 5 is complex (cf. 
example problem 5 d), then It may not appear possible to satisfy (76). (in 
problems where this occurs, the matching rule in the conventional method of 
matched asymptotic expansions also does not work unless one resorts to finding 
an intermediate expansion (cf. ref. l) . In the present method, although the 
relation (72) is valid and could be used if 5 were known completely, the 
relation (76) in that precise form may not work if 5 is larger than a 1 
because significant terms do not appear in f no that must appear in the 
asymptotic limit of f nx , and f nx can then not match asymptotically to 
(f no ) nx . The correct order of magnitude of 5 , and consequently the proper 
form to replace f no in (7 6), can then be determined by using the concept of 
the intermediate expansion. A displaced intermediate -expansion solution 
replaces f no in the matching rule and is determined by setting y = e 0 ^ - 5 , 
f° = ePf; specifying p in terms of a so that w is o(l) in the inter- 
mediate region; and determining a so that = e a (where 6 = cr^A) an( ^ so 

as to make possible the matching of the intermediate displacement expansion 
with the inner solution. Since y = 0 (b) in the "overlap” region where the 
matching is to be accomplished and for which an intermediate expansion is 
needed, this procedure makes A = 0 (l) in "y-space" in the intermediate prob- 
lem and determines the y-space where the overlap is possible. Then the 
large complex displacement, which caused the need for the intermediate expan- 
sion for matching, is determined in the matching process. In example prob- 
lem 5d, where a 1=6 and w = df/dy, it is found that p = a = 1/2. 

In classes of problems where the approach of the individual terms of the 
inner solution to the respective terms of the outer solution is not necessar- 
ily exponential, the asymptotic relation (72) can be replaced only by 

f nx ~ (f no ) ni + smaller alg as Y -> 00 ( 77 ) 

where Y = 0 is the locus of the nonuniformity and where "smaller alg as 
Y -► °°" means algebraically decreasing terms that are smaller as Y -> 00 than 
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those contained in (f no ) n:L . The appearance of the algebraically decaying 
terms in ( 77 ) (e.g., for class l) , as opposed to the exponential approach in 
class i s caused by the singularity at y = 0 in the outer expansion. The 
negative powers of y in the outer expansion, when transformed to inner vari- 
ables, introduce terms of lower order in e_ in the terms to be matched by the 
terms of the inner expansion. Each successive term of higher order in e in 
the outer expansion contributes terms to the lowest and succeeding orders in 
the inner solution. In example problem 1, where oi = l/o^ = e c = e 1 / 2 , the 
first term of the inner expansion is 

f ii = 6 -i/ 2 [-_y + (y 2 + 2 ) 1/2 ] = e“ 1 / 2 [Y -;L - (2Y 3 ) -1 + (2Y 5 )" 1 + . . .] 

which can match to (f 10 ) 11 = e -1/2 Y -1 only to within algebraically decaying 
terms as Y 00. Only additional higher orders in the outer expansion of f 
can account for these additional terms in the expansion of f 11 for large Y. 
The matching rule (77) allows for the extra algebraic terms. 

A form of the matching principle, which is analogous to Kaplun* s matching 
rule (45) in the conventional inner and outer expansion method, can be used in 
the present method in any of the classes of problems discussed above. If 
(f ni ) no (defined after eq. (74)) is written in inner variables and then trun- 
cated at n terms, one has [ (fni)no ]ni # Then a result corresponding to 
equation (45) is 


j*^ni^no]ni _ ^no^ni (yg) 

For problems in which 6 = 0 (including all of class l), equation (78) is 
equivalent to (45)* As in equation (45) > the displacements occurring in prob- 
lems of class 5 are not retained in the matching by equation (78). The form 
(78) is conveniently used to evaluate constants of integration in a closed- 
form inner solution, but if explicit asymptotic outer boundary conditions for 
the inner problem must be known before the differential equations in the inner 
problem are solved, they are conveniently provided by the asymptotic form of 
the matching principle given by (76) or (77)* It was noted above that if f° 
contains more than one displacement in class 5j then (76) does not apply. In 


that case, either (45) or (78) is suitable. 

The important results of this section are summarized as follows: 

In classes 1, 2, 3: f ni = (f°) ni (79a) 

In classes b, 5: f nl ~ (f°) ni + exp as Y -*• « (79b) 

In classes 1, 2, 3, f n ^ ~ (f no ) n i + smaller alg as Y -*■ 00 (80a) 

In class 5 with conditions (75) : f 11 '*' ~ (f no ) n ^ + exp as Y -* 00 (80b) 

In classes 1, 2, 3, 4, 5: [(f ni ) n °] ni = (f no ) nl (8l) 


Worthy of special note is the fact that, even for a problem in which it 
can be shown that a boundary -layer solution must go exponentially to an outer 
solution ((791) is satisfied), the respective terms may approach algebraically 
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to the terms of the outer solution according to (80a) if the outer expansion 
contains negative powers of y. This occurs in class 4 , of which two examples 
are provided above (p. 9)* 


The Composite Solution 

The outer and inner solutions may be combined, if desired, into a single 
composite expansion that is uniformly valid to appropriate orders in both the 
outer and inner regions, A rule for forming the composite solution is 
directly analogous to equation ( 46 ). It may be stated as 

fnc - pni + fiio _ (fno)ni (82) 


where . f nc is the composite expansion to nth order. Note that f no and 
(f no ) n=L in this method have a different meaning than in equation ( 46 ). 


Matching of Several Functions 

In a given problem one may wish (or need) to expand more than one 
function in asymptotic series and to match the several inner and outer func- 
tions. A singular perturbation problem may, of course, consist of a system of 
equations for which the order is reduced in the outer limiting process. Even 
in a problem with only one equation for one dependent variable f, the deriva- 
tives f y> f yy> etc., may be considered as different functions. Each such 

function, say g, should match according to the matching principle, with f, 
f°, F, etc., replaced by g, g°, G, etc., in the above development, and with 
era replaced by a g . When g is some function other than a derivative of f, 
0 g can be determined from the governing equations. 


Extension of the Method to a Class of Problems With Outer and 
Inner Functions Nonanalytic at e = 0 

If it is suspected or found that nonpowers of € c will occur in some 
higher orders in the procedure outlined above, then the inner and outer func- 
tions can be redefined as the sums of convergent series like f* in equa- 
tion ( 49 )* For the purpose of matching by (76) in class 5, the class of 
problems considered should then be restricted to those for which the first N 
terms in the expansions, which are to be determined, contain only powers of 
€ c . For the (N + l)st and higher order terms, (76) may not be applicable. 

Consider a continuous function f(x,y,e) with the asymptotic expansion 
( 47 )^ where (p ± = 1 and cp n (e) = o[cp n ^ 1 (e)] as e 0 and with the additional 
requirement on the asymptotic sequence {cp n } that for each cp n (e) there is some 
m > 0 such that 

€ m = o[cp n (e)] as € -> 0 (83) 

Thus, the sequence {cp n } may contain terms with factors like e 3 (log e)’” 1 or 
e 6 log e, but excludes terms like, for example, e~ 1/G . As in equation ( 49 ) * 
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let f* "be the sum of the series of terms cp n f n only: 

00 

f * = X c Pn( £ ) f n( x ^y) ln D *( x >y) ( 8i 0 

n=i 

which converges for sufficiently small e in some domain D*(x,y). (it can 
he shown (eVg . , see ref. 32) that if (47) converges at some (x,y), such as 
where a boundary condition is applied that either has a finite number of terms 
or is a convergent inf inite -series expansion of the same form as in (47 ) , then 
(47) converges in some neighborhood D*(x,y) of that location.) Then 

f( x ,y,e) = f*(x,y,e) + q.*(x,y,e) (85) 

where 

q *(x,y, e) ~ 0 + exp as e -> 0 (86) 

A convenient illustration of a typical function f* can he made as in 
equations (51) to (53) by adding a term e 4 log e/(l + y) to the right side of 
equations (51a) and (51c). In terms of y = eY, equations (52a) and (53b) 
would then have an additional term log e/(l + eY) on the right side and 

equation (52b) would have an additional § ( -l) ^Y J (e 1 ^ 4 log e) on the right. 

For the asymptotic expansion method, outer and magnified inner functions 
that are sums of convergent series in some domain are used (like f* and F in 
equations (51) to (53) with the additions noted above). Thus, it is assumed 
that 



and that 
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where 


f(x,y,e) 


Y 

F(x,Y,e) 


®i = 1 , ®n( £ ) 

q 1 (x,Y J ,e) 


f i (x,Y,e) + q. 1 (x,Y,e) 

-\ 

yAi(e) 


f 1 (x J ,Y J e)/a 2 (e) 


00 
r 1 


> f n (e)F n (x^Y) in some 

D i (x,Y) 

n=i 


°C a> n-i( € )^ e m = o[$ n (e) ] , 

m > 0 as € -> 0; 

0 + exp as e -* 0 with 

x,Y fixed J 


( 88 ) 


The remainder of the development is the same as in the preceding subsections 
(pp. 1 6 ff) except that the asymptotic sequences {cp n (€)] and [$ n (e)} must be 
determined by some means as in the conventional inner- and outer -expans ion 
approach (ref. l) . Since the satisfaction by f° of a boundary condition on 
f determines the displacement 5(x,e) (see subsections entitled "Definitions 
of inner and outer functions/ 1 pp. 21-25; and "Displacement Variables," 
pp. 27-29), the displacement may have the expansion (cf. eq. (64)). 


5 = cr s (e)A(x,e) 


a 5 (e) 


00 

9 n ( £ ) A n( x ) 

n=i 


The displacement variables are then 


CO > 

y = y - 8 = cr x Y - cr & cp n (e)An(x) 


f°(x,jr,e) = f°(x,y,e) 

and the expansion of f° is 

00 

?°(x,y,e) = ^ 9 n (e)f n (x,y) 

n=i 


(89) 


(90a) 


(90b) 


(91) 
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From equations (87) and (88), equation (70) still applies and (q 0 ) 1 - q a 
contains only terms that are changed from order exponentially small as e -*■ 0 
with x,y fixed to the order of some $ n (e) as e -* 0 with x,Y fixed. 
These terms then decay exponentially to zero as Y -*■ 00 . Thus equation (71) 
still applies, and to order cr 2 (e )$ n (e ) in inner variables, equation (72) 
applies. Now let f no denote n terms of f°(x,y,e): 


and 


n 

fn ° = X (92a) 

m=i 


n 

f ni = cr 2 (e)^ $ m (e)F m (x,Y) 
m=x 



and let (f no ) ni denote the terms to order cJ2(e)® n ( e ) in equation (92a) after 
it has heen transformed by equation ( 90 a) to x,Y variables. Then the match- 
ing rule (76) applies if f no contains all the parts of f° that would 
appear in (fO)ni. that is, f° can be replaced by f no in equation (72) if 
this condition is satisfied. At least to some order N for which 


cp n (e) = $ n (e) = [e c (e)] n 1 for a11 n < N ( 93 ) 

this condition is satisfied if conditions ( 75 ) are met. If equation ( 93 ) is 

true for some N hut conditions ( 75 ) are not met, then either an intermediate 

expansion is needed to replace f no in equation (76) or else equation (77) 
applies (for n <N), as discussed under The Matching Principle. 

If for some N it is found that for n > N terms with nonpowers of e c 
occur, then at least Kaplun’s standard matching rule (eq. ( 45 ) or (78)) can 
be used. 


APPLICATION OF THE METHOD WITH DISPLACEMENT VARIABLES 
TO HYPERSONIC VISCOUS FLOW OVER A SPHERE 


The purpose of this section is to apply the method developed above to the 
problem of hypersonic viscous flow over a sphere at high Reynolds number. 

This application is limited to cases for which external (shock-generated) vor- 
ticity has only second -order (and higher) effects on the solution for the 
viscous flow near the body. It will he seen that this limitation requires 
that the shock density ratio he not too small. (A list of references on the 
viscous blunt -body problem is given in reference 30- ) 

It has heen pointed out (ref. 30 ) that the second-order effects of 
external vorticity and displacement are "kindred effects" and that they are 
"global in nature" and, as such, are much more difficult to calculate than 
other second -order effects. In the point of view of the present approach, the 
presence of external vorticity is responsible for higher order displacements , 
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so that the inner boundary for the outer problem is affected to each 
successive order of the small parameter. Using the present approach, the 
displacement effects are determined in the matching process and subsequent 
solution of the inner problems. 

After definition of the mathematical model for the problem to be solved, 
and subsequent application of the procedure to obtain the asymptotic expan- 
sions, the results will be compared with a numerical solution of the model 
prob lem. 


The Mathematical Model for Hypersonic Viscous 
Flow Over a Sphere 

This approximate application to viscous flow is not intended as a 
thorough treatment of the hypersonic blunt -body problem for all regimes of 
flight (see, e.g., analyses in refs. 27 and 36-46). The object here is to 
isolate the role of the matching principle as much as possible and, in partic- 
ular, the role of displacement in the matching. The mathematical model will 
therefore be kept as simple as possible so that the basic principles will not 
be obscured by too much detail. To achieve this, a number of important 
aspects of some regimes of flight in the physical problem are ignored in favor 
of the desired simplicity. The model must, however, retain the essential 
character that we are studying. Simplifications that do not change the essen- 
tial character of the problem are the assumptions of: (l) constant density 
and viscosity behind the shock wave, ( 2 ) local similarity near the stagnation 
streamline (includes shock wave concentric with the body), and ( 3 ) discontin- 
uous properties across the shock wave, determined by Rankine -Hugoniot equa- 
tions. Also, the low-density effects of velocity slip and temperature jump 
are specifically neglected (ref. 47). All these assumptions and approxima- 
tions have been used and discussed by various authors. Included in the char- 
acteristics of the problem that must be retained are: curved body, curved 

detached shock wave, and consideration of the full governing equations 
(Navier -Stokes) . In the present treatment, curvature effects will be retained 
implicitly (by a transformation of coordinates), and need not be considered as 
second order or otherwise. The gas immediately behind the shock wave can be 
assumed to be in thermodynamic equilibrium for evaluation of conditions there. 
The flow is assumed to be steady. 

The Navier-Stokes equations for conservation of mass and momentum are: 

v • (pv) = 0 (94) 

p(v • v)V + VP = mv 2 v (95) 

where V is the vector operator "grad, " p the local mass density, p the 
pressure, V the velocity vector, and \i the viscosity coefficient (assumed 
constant), and where the equation of state for the flow between the shock wave 
and the body is taken to be 

p = constant = p 2 ( 96 ) 
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(subscript 
symmetry) . 
the others 
conditions 


2 denotes a value immediately behind the stock on the axis of 
Because of equation ( 96 ), the energy equation is uncoupled from 
and is not needed in the present analysis (except in finding 
across the shock wave) . 


Consider axisymmetric spherical coordinates (r,cp) with r measured from 
the sphere center and cp measured from the forward axial streamline (fig. 3) . 

Let the radius of the spherical 
body be rq, and that of the 
detached concentric shock wave 
(near the axis) be r s = rq^l + d) . 
Denote free -stream conditions by 
subscript °°. Let the velocity 
components, respectively tangen- 
tial and normal to the body, be 
VoqU and V^v, so that u a nd v are 
dimensionless and V = \l u 2 + v 2 

is the velocity magnitude. Define 
a dimensionless pressure function 

P = (p - Poo)/poo v co 2 ; and denote the 
shock density ratio by k = poo/p^. 
The boundary conditions for the 
flow near the axis (taken to be 
the no -slip surface conditions and 
the Rankine-Hugoniot shock condi- 
tions) in terms of these dimension- 
less variables are: 



Figure 3.- Configuration for sphere in hypersonic flov. 

At 


r = r. 


V 


At r = r s = r, (l + d) : 


u = v = 0 


v = -k cos cp 
u = sin cp 
p = (l - k)cos 2 cp 


( 97 ) 

(98a) 

(98b) 

(98c) 


For complete determination of conditions behind the shock, one needs to 
include also in the Rankine-Hugoniot conditions the energy relation across the 
shock: 


At r = r s : 

h + \ V 2 = K, + | V oo 2 ( 98 d) 

where h is the specific enthalpy, and an equation of state 

h=h(p,p) (98e) 

that may be represented, for example, by equilibrium air charts or tables. 
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If equations ( 9 I+) and (95) are put in terms of the coordinates r and cp, 
conservation of mass (eq. (94)) is satisfied by definition of a dimensionless 
Stokes stream function ^ by 


r 



sin cp 


Hr = - v (^) 6in * 

For convenience, we define a transformation 

£ = log(r/r b ) 


(99) 


(lOOa) 


so that 

r KA 

r Sr " 


(100b) 


(Note that in terms of physical distance from the body surface, zvy, the new 
variable £ is 


£ = log(l + y) = y + 0(y 2 ) (lOOc) 

and so is the same as y to first order for small distances. Thus, the 
variable t, is convenient both as a spherical coordinate in the outer problem, 
since it will simplify the equations, and as a boundary-layer variable.) With 
a Reynolds number defined as 


R = 


P 2 V oo r b 


Mv 


1 / PcoVbN 

k ^ V Poo ) 


( 101 ) 


the components of the momentum equation (95) become 
uu^ + vu^ + uv - R -1 e”^[u^ 4- u^ + 2v C p - u(sin cp) “ 2 

+ (cot cp) Uq) + u^ = -kp^ (l02a) 

O -1 r 

uv cp + ■ u ’ R e + V^ - 2v + Vqxp - 2u cp 

+ (cot cp)(v C p - 2u) ] = -kpj, (l02b ) 

where subscripts cp and £ denote partial differentiation with respect to the 
corresponding independent variables. Equations (102) may then be combined by 
cross -differentiating and by equating p^ to p ^ to eliminate pressure. 

Equations (99) with the transformation (100) may then be inserted to obtain a 
nonlinear partial differential equation for the single dependent variable \| / 
in terms of the independent variables cp,£ and the parameter R. The result 
is 
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tf^cp - 2 (cot cp)t^^ - + 2 (cot cp>|^ 2 + ^cpcpcp 

- 3 (cot cp)^^ + 3 (cot 2 cp - 1)^ - tcp*^ - 
+ + %^cpcp + (cot cpJt^ - 4 (cot cp)}^ 2 

= R -1 e^ sin cp[\[r^^ - + llf^ - - 2 (cot 

+ 2 %cpcp " 6t £qxp + 6 ( cot ^cp " (9 cot cp + 3 cot 3 cp)^ 

+ (8 + 3 cot 2 9 )f w - 2(cotcp)t q3qxp + ^ wq) ] (103) 

The boundary conditions for \|r, obtained from equations (97) and (98), are: 

At £ = 0: 

t; = ''f'cp = 0 (104a) 

At i = 5 S = log(l + d): 

i|f = k(l + d) 2 sin 2 cp cot cp (105a) 

\|/£ = (l + d) 2 sin 2 cp ( 105 b) 

and 


- (cot 9)^ 2 - + 'lyl^ “ vR -1 (l + d)sin cpl>^ 

- 3t^ + 2t ? + t 5qxp - (cot cp)t ?cp -2^ + 2 (cot cp^j 

= 2k(l - k)(l + d) 4 sin 4 cp cot cp (l 05 c) 

where one can take either v = 0 or v = 1 for the model (discussed he low) . 
Without loss of generality one can, of course, replace the condition a)t = 0 
in conditions (104a) with 

£ = 0 , f = 0 (104b) 

Condition (l05c) with v = 1 was obtained by first differentiating 
condition (98c) with respect to cp and combining the result with equa- 
tion (102a). Thus, v = 1 implies strict application of the Rankine-Hugoniot 
shock conditions (eq. (98)). However, since those conditions are derived 
assuming inviscid flow on both sides of the shock, one might just as well take 
v = 0 (ref. 48). The condition with v = 1 is more consistent with the 
Rankine-Hugoniot conditions, but the condition with v = 0 is simpler. In 
any case, the value of v does not affect the perturbation problem in the 
first- and second-order solutions (refs. 46 and 48). (Some subsequent work 



shows that this is time only if, as assumed here, the vorticity interaction is 
not greater than a second -order effect.) Moreover, the differences in results 
calculated for the two values of v in the exact solution could he taken as 
an indicator of the influence of neglecting the shock-wave -viscosity effects 
on the flow behind the shock. 

Equation (103) and boundary conditions (104) and (105) contain two 
independent parameters, k and R. The dimensionless shock detachment distance, 
d, is an unknown function of k and R, to be determined in the solution. 


Solution of the Problem by Assumed Power-Series 
Expansions With Displacement Variables 

The problem of equations ( 103 )* (104), and (105) is to be solved 
asymptotically for \|r(cp,£,k,l/R) as l/R -> 0. The parameter k is assumed to 
be nonvanishing as l/R -> 0. The result of substituting 

i|r(<p,£,k,l/k) ~ t x (cp^A) as l/R -* 0 ( 106 ) 

into equations (103) > (l04), and ( 105 ) is that the order of the partial 
differential equation is reduced and the no -slip condition, \|r^ = 0 on £ = 0, 
cannot be satisfied. This problem is known (and will be seen) to belong to 
class 5, as defined above (pp. 5-9) , and will be treated as a problem of that 
class . 

The outer function representing the stream function . - According to the 
method outlined above (cf. eq. (55a)), let 

k,i/R) = t°(cp,Ok,i/R) + q°(cp,OkA/R) (107) 


and assume 


d(k,l/R) = d°(k,l/R) + q d (k,l/R) ( 108 ) 

where i|f° and d° are assumed to he analytic with respect to some e = e(l/R) 
at l/R = 0 so that 

t°(cp,£,k,l/R) = f 1 (9,Ok) + et 2 (<p,Ok) + e 2 t 3 (cp,Ok) + • • • (109) 

d°(k,l/R) = di(k) + ed 2 (k) + e 2 d 3 (k) + . . . (110) 

and where 


q° ~ 0 + exp 
q d - 0 + exp 


as l/R 0 with 

as l/R 0 with 


cp,£,k fixed 
k fixed 


(in) 


Since the limit function in ( 106 ) was found to he a nonuniform represen- 
tation of at £ = 0, hy equation (l 07 )> the outer function i|r° represents 
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\|r to within exponentially small amounts (q°) as e -*■ 0 everywhere except at 
i = 0. Then substitution of equations (107) and (108) into (103), (104) , and 
(l05) shows that i|r° must satisfy equations (l03) , (104), and (105) with d 
replaced by d° and with t replaced by \|/° at all £ except very near 
£ = 0. According to the method outlined above, then, since wi = is the 

function that cannot satisfy a boundary condition at £ = 0, that condition is 
not required of w° = but the condition \|r = 0 can also be required of 

\|r° if it is allowed to b be satisfied at a displaced boundary £ = 8 (see 
p. 23), with the displacement 6 to be determined by compatibility with the 
inner solution. Assuming that 5 is real in this problem (see subsection 
"Displacement Variables," pp. 27-29), we can write 

5/ffi = A(cp,k,l/R) = Ai(cp,k) + eA 2 (cp,k) + e 2 As(cp,k) + . . . (112) 

where c?i = 0i(l/R) (to be determined) represents the order of magnitude of 
the width of the nonuniform region measured in units of £ . The problem for 
\|/° has now been derived as: equation (103) with \|r replaced by i|r°; and the 


boundary conditions: 


5 = 01 A , 1r° = 0 

(H3a) 

i = log(l + d°) : 


= k(l + d°) 2 sin 2 cp cot cp 

(1131) 


i|f£° = (l + d°) 2 sin 2 cp 

(113c) 


and 


[condition (l05c) with i|r replaced by i|r° and d by d°} (ll3d) 

The quantities e = e(l/R), 01 = 0i(l/R), A = A(cp,k,l/R), and d°(k,l/R) are to 
be determined. Note that 01A is a displacement of the independent variable 
£ and is therefore proportional to the actual displacement distance to first 
order only (see eq. (100c)). 

A dimensionless magnified viscous-flow displacement. A*, of the y 
coordinate may be defined, with use of equation (lOOc), by 

A' = y s /0i = (e 5 - l)/0i (114) 

The displacement distance is related to the shock standoff distance as follows. 
Denote by = r^(y s + l) the radius of the frictionless body that would 
produce the flow field that corresponds to the actual outer flow when viscos- 
ity is present (cf. ref. 33)- The shock standoff distance from the equivalent 
frictionless body is denoted by r^d* (fig. 3); from the actual body, it is 
r^d. (For v = 0, d* = di (given by the inviscid solution of Lighthill) . 

For v = 1, d* will depend on R, but not to the order to be considered here, 
that is, d* = di + v0(l/R) . However, the complete solution for d* is shown 
at the end of this section for the comparison of the asymptotic expansion for 
A* with the exact solution of the model problem.) The shock radius is then 
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r g = rgl + d) = r 5 (l + d*) 


(115) 


The dimensionless magnified displacement is then given by 


A* 


^5 = r S I r b = d - d* 

(7i air b ohl 1 + d*) 


( 116 ) 


Equating (ll4) and (ll6) gives the relationship between d-d* and 5: 

d - d* = (l + d*)(e 6 - 1) ( 117 ) 

Note now from equation ( 117 ) that, since d* and 5 are analytic with respect 
to e c at e = 0, so also is d. Then q^ = 0 in equation (108) and 

d° = d = di(k) + eds(k) + e 2 d 3 (k) + . . . (ll8) 


The assumptions of the concentric shock and local similarity of the flow 
solution require the dimensionless standoff distance, d, and r^ both to be 
functions only of k and R, not varying with cp. 

The displacement variables cp and \ are defined (cf. eqs. (64) and ( 67 )) 
by: 

$ S 9 , C = 5 - ffi A (119a) 

The outer displacement dependent variable (cf. eqs. ( 68 ) and ( 69 )) is defined 

by 

t°($l,k,l/R) = i|r°(cp,£,k,l/R) (H9b) 


= + e\f 2 (cpl,k) + e 2 t 3 ($,Lk) + . . . (119c) 

Substitution of equations (109)* ( 112 ) , and (ll8) into the differential equa- 
tion for i|f°, ( 103 ) , and into the conditions ( 113 ) and letting l/R -> 0 (so 
that also the as yet unknown quantities e and a x vanish) gives the problem 
for \|/ x * In finding the boundary conditions for "vjr x it is assumed that 
is analytic with respect to £ at £ = CiA and at £ . The condition at 
£ = o*iA is then transferred to £ = 0 by Maclaurin* s series, from which 

t°(cp^criA) = \|r°(cp,0) + cr 1 At^ o (cp,0) + (l/2! ) (<JiA) 2 t^(cp,0) + . . . = 0 (120) 

and into which equations ( 109 ) and (112) must be substituted. Similarly, the 
conditions at the shock [£ = £ s = log(l + d) ] are expressed in terms of 
conditions at 


(121a) 
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£ S1 = log(l + da.) 



where 


L - 


5 fll - e 


1 + d x 


+ e* 


1 + dx 2 VI + dx 
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by means of Taylor’s series. 

Solution of t he first oute r problem . ~ The resulting problem for 


♦ l5 + l55q , - 2(c°t - * l5 + l5q) + 2 (cot cp>lf^ 2 + t 1 ^ 1 cpcpcp 


- 3 (cot cp)t 1 ^ lC p C p + 3 (cot 2 cp - 1)^^ - - ^ iqX p^ i<p 

+ 5+^^ + ^ iCp t lW + (cot <p)t iq5 t M - Meet cp)t icp 2 = 0 
with the conditions 

ti(<P>0) = 0 

^(<P^ S1 ) = (l + di) 2 sin 2 cp 
't r 1 cp( c P^sx) = k ( 1 + d i) 2 sin - 2 <P cot 9 

and 


~ (cotcp) ^ 2 - +x9 + i55 + VisL, 

^ ^S1 

= 2k(l - k)(l + d x ) 4 sin 4 cp cot cp 
This problem has the separable solution 

= (sin 2 9)fi(0 

where fi(£) is found from the separated ordinary differential equation 

fi m ~ 5f i" + 2f i* + 8f i = 0 

(where ( ) 1 generally denotes differentiation of a function with respect 
its argument) and has the form 


(121b) 

is 

(122a) 

( 122b ) 
( 122c ) 
(l22d) 

( 122e ) 

( 123 ) 

( 124 ) 
to 
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(125a) 


f = Cne 4 ^ + Ci 2 e 2 ^ + Ci 3 e ^ 

= C u (r/r b ) 4 + Ci 2 (r/r b ) 2 + Ci 3 (r/r b ) 1 ( 1251 ) 

= aio + an£ + a 12 ^ 2 + ai 3 £ 3 + . . . (l25c) 

This is the well-known solution of Lighthill (ref. 49 ) for the inviscid flow 
over a sphere. The boundary conditions determine that 

Cxi. = (1 - k) 2 /lOk(l + dx) 2 

Ci 2 = (4k - l) /6k / (126a) 

C 13 = (1 - k)(l - 6k) (1 + di) 3 /l 5 k j 

and that the function d x = di(k) is the solution of the algebraic equation 
obtained from 

Cn + + C 13 = 0 (126b) 

The constants in equation (l 25 c) will be needed for matching to the inner 
solution and are easily found to be: 


aio = 0 (127a) 

an = 4 Cn + 2 C 12 - C 13 = ^ ([ ' Vdi ' P + 4 k - 1 (l 2 Tb) 

a 12 = 8Cn + 2C 12 + ^ C 13 = + “ 1 (l27c) 

- f C X1 + ^ 0l2 - \ C 13 - ^ + 3 ( 4 k - 1 )_ 


a n a i2 

* ' T + 5 T 

The significance of the constants an and am can be observed from the fact 
that the dimensionless velocity u and the dimensionless vorticity 
w = (r- t) /V 00 )(y X "\f) at the body surface in the first-order inviscid solution 
are : 

u ll3 = &n sin cp , w lb = ( a n “ 2a 12 )sin cp (128) 



45 



To find the terms for the higher order outer solutions etc.) the 

parameters e = e(l/R) and cri = 0i(l/R) must he known. Therefore, the inner 
problem (of determining the inner function \|A, to be defined) must be formu- 
lated. Outer boundary conditions for the inner problem will need to be 
supplied by the matching principle which, in the present approach, uses the 
displacement variables. For that purpose, the first term in the outer dis- 
placement expansion (ll9c) is found as follows: From equations (l09)> (ll9b) , 

and (ll9c). 


= \|/ 1 (cp,^,k) + 0(e) (129) 

Into the right side of equation (129) we then substitute (123), (l25c), and 
the transformation (l!9a) to obtain 

ti = sin 2 cp(an£ + a 12 £ 2 + a 13 l 3 + . . .) (130) 

The inner function representing the_stream function.- The inner function 
ijf 1 , which represents ^ , at least asymptotically, in the inner region as 
l/R 0, is found hy magnifying the variables in such a way that the nonuni- 
formity of the problem (reduced-order differential equation with loss of a 
boundary condition as l/R -> 0) is removed. Following the method outlined 
above, one can define magnified variables Z and if such that (cf. eqs. (60 )) 

Z = 5/cri (131a) 

and 

?(cp,Z,k,l/R) = (l/a 2 )^/ i (cp,Z,k,l/R) (131b) 

where 

i|r(q>,5,k,l/fe) = t 1 (cp,Z,k,l/R) + q i (cp,Z,k,l/R) (l32) 

and ¥ is assumed to be analytic with respect to some € = e(l/R) at l/R = 0 
so that 

l(cp,Z,k,l/R) = ^i(cp,Z,k) + e¥ 2 (<p,Z,k) + e 2 ¥ 3 (cp,Z,k) + • • • (133) 

and where 

q-'-(cp,Z,k,l/R) ~ 0 + exp as l/R -»• 0 with <p,Z,k fixed (134) 

The quantities e = e(l/R), 01 = 0i(l/R), and a 2 = a 2 (l/R) are to be 
determined. 

Substitution of equations (l3l) through (134) into the boundary 
conditions (l04a) and (l04b) shows that H must satisfy the conditions 

T=¥ z =0 at Z = 0 (135) 

To determine oi and cr 2 , consider the rules stated under Magnifying 
Factors. First, since the function that does not satisfy one of the inner 
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■boundary conditions in the outer solution is w° = ^°, we require that 

he o(l) in the inner region (as l/R -► 0 with Z fixed). Then, from 
equations ( 131 ) 

= (cr 2 /ai)^ z = 0(l) as l/R -* 0 with Z fixed 

and, since Y% = 0(l) as l/R -> 0 in the inner region, we take, without loss of 
generality, 

0*1 = 0*2 = o' ( 136 a) 

If equations ( 13 I) through (134) are substituted into the differential equa- 
tion (103) 9 and if the equation is then multiplied by a, and l/R is made to 
go to zero, application of the principle of least degeneracy requires that 

o^R = 0(1) as l/R 0 


so that 


cr = R _1/2 (136b) 

The inner problem is the problem for Y(cp,Z,k,l/R) which is composed of 
conditions ( 135 ), additional outer conditions to be supplied by the matching 
principle, and the differential equation (l03) with equations ( 131 ) and (l 32 ) 
substituted. Except for the outer boundary conditions, the inner problem for 
Y contains only integral powers of l/R. However, the outer problem contains 
cr = R -1/2 in its inner boundary condition, and it will be seen that R“ 1/2 
also occurs explicitly in the outer conditions supplied to the inner problem 
by the matching principle. Hence, as described above, we then assume that 

e = R” 1/2 (137) 

is the form of the small parameter with respect to which \| r° and Y are 
analytic at e = 0. 

Solution of the second outer problem . - The problems for the respective 
terms of the outer expansion (109) are determined by substituting equa- 
tions (109) j (llO), and (ll8) into the problem for \|r° as stated in the sub- 
section "The outer function representing the stream function, ” pp. 4l-4^. The 
problem for ijr . l is given by equations (l22) with the solution given by equa- 
tions (123) through (126). The problem for \|r 2 contains a great deal more 
terms than the problem for \| r and so will not be written out here. Note, 
however, that the condition on \|r 2 from equation (ll3a) is 

t 2 (<P>0) = -A i ^ 1 ^(cp,0) 

which contains the unknown A 1# The problem for can be solved at this 

point if the form of is known. To continue the assumption of local 

similarity in the higher-order outer solutions, we need to have the displace- 
ment surface concentric with the body near the axis (i.e., A = constant), 
like the shock wave. Let us carry on with this assumption to find the second 
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outer solution,, but then solve the first inner problem, without making the 
ass ump tion a priori there, to determine the actual significance of assuming a 
concentric displacement surface (i.e., of taking Ax to be independent of cp 
near cp = 0). The second outer problem has the separable solution (found in a 
different form by Van Dyke, ref. 30) : 

t 2 = (sin 2 cp)f 2 (0 (138) 

where f 2 (£) is the solution to the ordinary differential equation 

f|* - 5f” + 2f^ + 8f 2 = 0 (139) 

Equation (139) has the same form as (124) , and so has the solution 

f 2 (C) = C 2 xe 4 ^ + C 22 e 2 ^ + C 23 e'^ (l40a) 

= a 2 o + 3-2 1C t a 22 £ 2 + . . . (l40b) 


The boundary conditions for and those for lead to determination of 

the constants C 21 , C 22 , C 23 , and d 2 . They are found to be 


* 

C21 = “(I - k) 2 d 2 /5k(l + dx) 3 , C 22 = 0 

C 23 = (l - k)(l - 6k)d 2 (l + dx) 2 /5k 


(lUa) 


and, -with 


C21 + ^22 + ^23 “ 


(l4lb) 


d 2 is determined explicitly as 


d 2 = (l + dx)Ax 


(l4lc) 


(Note that this result can be obtained directly from expanding equation (117) 
and using (112) and (ll8) . When A 2 is known, one similarly obtains the 
further result that d 3 = (l + dx)(A 2 + Ax 2 /2).) For a given k, the quantity 
dx in equations (l4l) is known from the first outer solution. The value of 
Ax must be obtained from the first inner solution. For use in the matching, 
the constants in equation (l40b) are found to be 


a 2 o - -axxAx > &2X = (2axx " 2ax2)Ax 

app = (axx “ 3a-x2)Ax > a 23 = (l/3)(7a-xx ~ 13a-x2)Ax 


(142) 


The same procedure that was used to obtain equation (130) from one term of 
\|r° now leads to two terms of V° corresponding to two terms of 7° • 

t° = sin 2 cp{(aix^ + ax2^ + aaaC + . . .) 

+ e[(2a 12 A 1 + a 2i )£ + (3^x3Ax + a 22 )£ 2 + . . .] + 0(e 2 )} (143) 
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Upon comparing equation (l43) with (ll9c) we see that t 2 is the factor 
multiplying e in (l43). 

Oute r boundary conditi ons for the inner problem . - Equation (l43) is used 
to find the outer "boundary conditions for the inner problem. The terms to 
order e in equation (l43) constitute i f 2 ° (cf. eq. (73))- From equa- 
tions (112), (119a), and (l31a) we have 

cp = cp, \ = e(z - Ai) - e 2 A 2 - € 3 A 3 "... (l44) 

which, when substituted into f 2 °, lead to 

(t 20 ) 21 = sin 2 cp{eaix(Z - A x ) + e 2 [ai 2 (z - Ax) 2 

+ (a 2 x + 2ax2^i)(Z “ ^i) ” axx& 2 ]} (l45) 

as defined in the discussion following equation (74) • The matching principle 
(76) requires that 


i|/ 2 i ~ (i}r 20 ) 2 l + exp as Z -> °° (l46) 

where 

^2i _ + € 2 y 2 (147) 

(cf. eq. (74)). From equations (l45), (146) , and (l47) we then find the outer 
boundary conditions for ¥x and ¥ 2 : 

^x ~ sin 2 cp[axx( z - A x ) ] + exp as Z -*• °° (l48a) 


Wq ~ sin 2 cp[axs(Z - Ax) 2 + (a 2 x + 2ax2^i)(Z 


(l48b) 


So luti on of the first- ana second -ord er inner problems .- Now substitution 
of equation (.133) Into the partial differential equation for l(cp,Z,k,l/R) 
found above (see discussion following eq. (136b)) and into the conditions 
(l35)> and letting e -> 0 successively, leads to the following problems for 
the first two terms in the inner expansion (133) : 


c>Z 


¥1 = Y lZ = ° at Z = 0 


*lZ*x<pZ - ViZZ ~ ( cot ^^xZ 2 " (5ln ^iZZZ 


= 0 


(149) 


with equation (l48a) as the outer condition, and 



_s_ 

dz 


-1^ 


Wapz + V^z ' W*zz - Wap - 2 < cot ^Z f 2 z - < sl " *2ZZZ_ 


= f _¥ _ - 2 (cotcp)'F „ z - 51 I __ - sin cp 

lZ lZcp Y lZ 19 xZZ Y 


^lZZZ Z hzZZZ 


> (150) 


¥ 2 = ¥ z =0 at Z = 0 


with equation (l 48 b) as the outer condition. The problem for ¥1 is, of 
course, the same as would be obtained by Prandtl 1 s first-order boundary-layer 
theory. 

As discussed above (prior to eq. (138)), in solving the problem for ¥1 
we will not assume a priori that Ai = constant (not varying with cp), but 
■will use the equations that determine the form of Ai and then observe 
the significance of requiring A± = constant. The bracketed quantity in 
equation (l^ 9 ) must be a function only of cp, say 0 (cp) • Differentiation of 
condition (l^a) to obtain asymptotic expressions for the various terms in 
cp(cp) then leads to 


Now assume 


where 


2 

§(cp) = an sin 4 cp cot cp 


n =-s/2F Fi(ti) 


pep 

I = l(q>) - an / sin 3 cp dcp 


- an 


1 - cos cp - ^ (l - cos 3 cp)~l 
■X J 


an sin 4 cp 
4 




1 + ir sin 2 cp + 0(sin 4 cp) I 

3 A J 


(15D 


(152) 


(153a) 


and where 


T) = T| (cp_y z) = AZ 


(153b) 
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(153c) 


, , / . an sin 2 cp 7 

A = A(sxn cp) s — : = 2a u 

V2| 


1 - ^ sin 2 cp + 0(sin 4 cp) 


Also define a function A by 


A = A(sin cp) = C ? S 4 C ^ 
an sin 4 cp 


1 

2 


(154) 


1 - g- sin 2 cp + 0(sin 4 cp) 

The problem for ¥1 given by equations (l48a) and (l49) then becomes 

Fi»tt + FiFi" = A[(Fi') a - 1] 

Fi(0) = Fi'(O) = 0 
F x ~ T] - AA X + exp as rj 00 


( 155 a) 

(1551) 

(155c) 


or 

Fi* - 1 + exp as r\ 00 ( 155 c 1 ) 

where ( ) f here denotes differentiation with respect to rj . A local 
similarity solution is obtained by assuming A to be a constant appropriate 
to the particular value of cp (cf . Stine and Wanlass, ref. 50j or Lees, 
ref. 51 )- In particular, for the region near the stagnation streamline, 
neglect of (l/6)sin 2 cp in comparison to unity in equation ( 154 ) leads to a 
locally similar solution with A = 1/2. Further, neglect of (l/6)sin 2 cp in 
comparison to unity in equat ion (l53c) then leads to a constant value for A x 
near cp = 0 with A = N /2an. For consistency, that is, in order to neglect 
all quantities of 0(sin 2 cp) as cp ^ 0, we can also neglect terms as small as 
(l/3)sin 2 cp in equation (l53a), so that we now take 


| = (1/4) an sin 4 cp ; A = N /2a 11 ; A = 1/2 ( 156 ) 

in the determination of ¥ x by equations ( 151 ), (l 53 "b), and ( 155 ). 

We see that the assumption of local spherical symmetry in the outer 
solution and further local similarity in the inner solution require neglect of 
0 (sin 2 cp) in comparison to unity in some quantities in the inner solution 
(cf. discussion by Kao, ref. 46, p. 1893) • This then yields constant A x 
which, further, will allow the condition of local spherical symmetry in the 
outer solution to continue to second order. 

Equation (155&) is equivalent to the equation derived by Falkner and 
Skan (ref. 52). The numerical solution corresponding to A = l/2 was first 
given by Homann in 1936 (see ref. 53)- More general solutions were studied 
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"by Hartree (ref. 54). A convenient tabulation of the solution for A = 1/2 is 

given in reference 55 (table V.3, p. 237) where , for example, 

' 

Fi"(0) = 0.9278 

and l (157) 

Pi = lim (r) - Fx) = 0.8046 

■^->00 

From this solution 

Ax = Px/A = 0.8046 /72aii (158) 

where an, depending on k, is given by equation (l27b) . This agrees with 
Van Dyke's evaluation (ref. 30) of the first - order displacement thickness 
(there referred to as the displacement thickness) at the stagnation point for 
the special case k = l/6 for which axi = 2/3- 

The second-order boundary -layer problem is given by equations (150) with 
the outer condition (lijSb) . A local similarity solution near cp = 0 is 
obtained as 

= (l/A) j~2£ F 2 (t]) (159) 

(with neglect of 0(sin 2 cp) in comparison to unity, as in the problem for ¥ 1 ) 
where | and A are given by equation (156) and t| = AZ. The resulting 


ordinary differential equation is 

(d/dr)) [F 2 '" + FxF 2 " - Fi'F 2 ‘ + Fx"F 2 ] = 5FxFx" + 6 Fi»" - t)Fx"" (l60a) 
with the conditions 

f 2 (o) = f 2 » (0) (160b) 

F 2 ~ (ax2/ a n) (7 - P1) 2 + 2Px(q - Pi) - p 2 + exp as 7 -» 00 (l60c) 

from which also 

F 2 l - (2a 12 /an) (-q - Pi) + 2pi + exp as r\ ™ (l60d) 

F 2 n ~ 2ai 2 /an + exp as y] -> 00 (l60e) 

and where 

P 2 = 2anA 2 (l6l) 
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The quantity is the second-order term for displacement in equation ( 112 ). 

With use of equation (155a) -with A = l/ 2 , equation (l 60 a) can he integrated 
to 


= (5/3)(F 1 F 1 l - r\) + (ll/3)Fi n - nFi”» + C x ( 162 ) 

where L is the linear differential operator defined by 

L ' F ^ 3 + Fl W ~ Fl ' m + Fl ") F2 (l63) 

and where C± is a constant of integration determined by the asymptotic 
boundary conditions to be 


Cl = -(l/3)Pi (164) 

The problem for F 2 is most easily solved by first splitting F 2 into 
two parts as follows. Let 

F 2 (t}) = P 2 i(tj) + F 2 "( 0 )F 22 (t)) ( 165 ) 

where F 21 (t}) is defined as the solution to the problem: 

!{F 21 } = (5/3)(PiPi* - ti) + (H/3)Fi" - tjPi'” + Ci (166a) 

Fai(O) = F 2 i'(0) = F 21 "(o) = 0 (l66b) 

Then the problem for F 22 (t}) becomes 

L{F 22 } = o (167a) 

F 22 (0) = F 22 »(o) = 0 ; F 22 "(0) = 1 (167b) 


It is easily found by asymptotic integration that equations (l66a) and (167a) 
have the asymptotic solutions 

F 2 i ~ K 2 (t] - Pi) 2 + Ki(t) - Pi) + K 3 + exp as ^ ^ 00 

, ' (168 
F 22 ~ K 4 (t) - Pi) + K 5 + exp as t) -» °° 

where the constants K 2 , K 3 , K 4 , and K5 are determined by the initial condi- 
tions (l66b) and (167b), and where K x = 2 p x . The constants can be obtained 
by numerical integration of equations (166) and (167) and asymptotic evalua- 
tion of the limits (found from (168)): 
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K 4 = (1/2)F 22 H 


K 2 = (1/2)F 2 i(“) , 


K 3 = lim [F 2 i(t]) - 2p x (ri - Pi) - (l/2)(r) - Pi) 2 F 2 "(t))] 

Tf* 00 

Ks = lim [F 22 (r|) - (1/2) (n - Pi) 2 F 22 (t}) ] 

T |->00 

Then it is found from conditions (l60c), (l6od), and (l60e) that 


(169) 


and 


F 2 "(0) = lim 

T) - * 00 


2ax 2 /an - F 2 i(ri) 

F 2 i’(Tl) 


a ig/ a ii - K 2 

k 4 


p 2 = -Ks - K s F 2 "(0) 


(170) 


Note from equations (167) that the solutions for F 2i and F 22 depend only oh 
prior solution of the first-order houndary-layer problem (for F x ), and so are 
independent of the external conditions. For an arbitrary k, numerical compu- 
tation of the constants K 2 , K 3 , K 4 , and K 5 then leads to: 


F 2 "(0) = 0.6790 56 (2ai 2 /a 11 ) - 0-774713 

Pg = 1.3796 - 0.701957(2ai 2 /an) 

■where a xx and a X 2 are given hy equations (127) and where d x 
equation (l26h). 


(171) 

is the root of 


Determination and matching of pressure . - As some difficulty in the past 
has heen associated with the determination of the second-order pressure, con- 
sider now the pressure and the matching of pressure. Corresponding to \|r° 
we have 


p°= p 0 (cp,£,k,,l/R) = p l (cp^,k) + ep 2 (cp,£,k) + e 2 p 3 (cp,£,k) + . . . 

(172a) 

where p° is the assumed analytic outer function representing p asymptoti- 
cally in the outer region (as l/R -> 0 with £ fixed). Corresponding to \jr°, 
we have 


P° = P° = P°($^Lk,l/R) 

= Pi($>Lk) + ep 2 ($,£,k) + e^C^Lk) + . . . (l72b) 

The inner function pi is to represent p asymptotically as l/R -> 0 with Z 
fixed. Thus, following the convention outlined above, define 
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P(<P>£AA/R) + exp 


(173) 




cr p (l/R)P(<p,Z,k,l/R) s p'*'(cp,Z,k,l/R) ~ 

as l/R 0 -with Z fixed 


where 


P = Pi(cp,Z,k) + eP 2 (cp,Z,k) + e 2 P 3 (cp,Z,k) + . . . (l r Jk ) 

is assumed analytic with respect to e at e = 0 . The factor a p (l/R) will 
he determined to make Pi the first nonzero term in equation (174). The 
inner and outer pressure functions will he seen to match according to 

p n i - (p 110 ) 11 ^ + exp as Z (l 75 ) 

as expected (cf. eq. (76)). 

For determination of the terms in equations (172) and (174), equa- 
tions ( 102 ) can he put in terms of i|r. The equations for p° are then the 
corresponding equations in terms of t| f° and the equations for p 1 are the 
corresponding equations in terms of 


The first- and second -order terms of the pressure function in the outer 
problem are then to he found from: 

k( 5 p 1 /Scp) = -sin cp cos cp e” 4 ^[(f! f ) 2 - 2fifi" + 2fif 1 1 ] (176a) 

-sin 2 cp - (f i* ) s ] - cos 2 cp e“ 4 ^[ 4 fif 1 ' - 8fi 2 ] (176b) 

and 


k tep 2 Ap) = 

k(dp 2 /S(;) = 


-2 sin cp cos cp e 4 ^[f 1 , f 2 * - fif 2 " - fi"f2 + (fifa)'] 

-2 sin 2 cp e“ 4 ^[f 1 f 2 * + f x 'f 2 - f x 'f 2 '] 

-4 cos 2 cp e 4 ^[(f 1 f 2 )' - 4 f x f 2 ] 


(176c) 


(I 76 d) 


Equations (176a) and (176c) may he integrated with respect to cp, then 
differentiated with respect to £ and equated, respectively, to equa- 
tions (176b) and (l76d) to evaluate functions of integration, making use of 
the differential equations ( 12 ^) and ( 139 ) • The constants of integration are 
evaluated hy the boundary conditions on fi and f 2 at (transferred to) £ S1 - 
The results are 


£ 

-kp 1 = sin 2 cp e -4S 


| (fi*) 2 - fifi" + fifi’ 


+ 2e -4 ^f i 2 + | k 2 - k 


( 177 a) 


55 



-kp 2 = sin 2 9 e~ 4 ^[fi'f 2 ' - fif 2 " - fi"f 2 + (fif 2 )'] + 4e‘ 4 ^f 1 f 2 (177*0 

In the equations for the inner problem it is found that, for nonzero Pi, 
cTp(l/R) must be 0(l) as l/R -> 0; thus, we take Op = 1. The equations for 
Pi are then 


k sin 2 cp[dPi/d( sin 9 ) ] = -$(cp)/cos cp = afi sin 3 cp 1 


kP^ = 0 


( 178 ) 


where ®(cp) is the bracketed quantity in equation (l49)j with the result 


kPi = -(l/ 2 ) afi sin 2 cp + constant 

Since k(p 1 °) li = -(l/2)a 2 i sin 2 cp + k - (l/2)k 2 and p 1 ^ = Pi, the matching 
rule ( 175 ) i £ satisfied for n = 1 by 

kPi = -( 1 / 2 ) afi sin 2 cp + k - (l/ 2 )k 2 (179) 

The equations for P 2 are (with neglect of 0(sin 2 cp) in comparison to unity 
and with use of equations (l 55 a) and (l 62 ) ) : 


k[SP 2 /ci(sin cp) ] = n/ 2 ail 3/2 (sin cp) [ (2/3) (Fi" + FiFi* ) + (l/3)i) + Ci] 

k( 8 P 2 / 8 r)) = (l// 2 )an 3 / 2 (sin 2 cp)(Fi ') 2 
which give 


(180) 


kP 2 = (l//2)an 3/2 (sin 2 q>) [( 2 / 3 ) (Fi" + PiPx* ) + (i/3)tj + C x ] + cor^sSt° (l8l) 
where Ci is given by equation (l64) . Since 

k(p 20 ) 21 = [-(l/2) afi sin 2 cp + k - (l/2)k 2 ] + e sin 2 cp[afiZ 


- ana 2 i + 2ai 2 a 2 o - ana 2 o] (l 82 ) 

and p 2i = Pi + eP 2 , it is found that equation ( 175 ) is satisfied exactly for 
n = 2 only if the constant of integration in equation (l8l) is zero. 

The body surface pressure is given by 


P ~ 5 *) = Ph ~ p 7 cp, 0 ,k, | J = Pi(<p, 0 ,k) + eP 2 (cp, 0 ,k) + . 


^ 00 ^ DO /1q 


\ 


= ( 1 ■ I ■ iir sln2 9 J + 


e(an) 


3/2 


42 k 


I Pi"(0) + Ci 


sin 2 cp + 


(183) 
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Results for velocity and skin friction .- The dimensionless velocity 
components u° and v° corresponding to the outer function \|/° are defined by 


u 


o — 


e" 2 ^£° 
sin cp 

. -2t,> r O 


- (sin cp)e 2 ^[fi* (£) + 6f 2 ’ (O + 0(e 2 )] 


v o = ~ e jg 

sin 9 


= -(2 cos cp)e 2 ^[fx(0 + ef 2 ({;) + 0(e 2 ) ] 


(l8 4a) 

(18 4b) 


The velocity components corresponding to the inner function ^ are 
-2eZ 


. e ¥ z 

u 1 = — : = sin cp 

sin cp Y 


auFi* (n) + e J~ 2 ^ ( f 2 '(ti) - 2t 1 F i '(r 1 )]+ 0(e 2 )l (l84c) 


-2€Z 
-e e% 


V 


£ 


sm cp 


= -e(cos cp) i N /2a^7 Fi(rj) + €[F 2 (t)) - 2rjF 1 (’n)] + 0(e 2 ) [ (l84d) 


It is easily shown 
and v (eq_s. (l84)) 


for n = 1 and n = 
Cf = (V^PoJoo 2 


that the inner and outer functions representing both 
satisfy the matching rules (cf. eq.. (76)): 

u ni „ (u no ) ni + ex p aS z -> 00 

v ni - (v 110 ) 11 ^ 4- exp as Z -> 00 

2. The skin-friction coefficient is 


= J^b _( 

' P»Vb V 



2 

kR sin cp 




u 


(185) 


d86a) 


2e 

k sin cp 


= e sin cp 


[¥ iZZ (cp,0,k) + eY 2ZZ (cp,0,k) + 0(e 2 ) ] 

^ 2aj -^ 3/2 *i”(0) + e % i F 2 "(0) + 0(e 2 ) 


= e sin cp [Cf + eCf 2 + 0 (e 2 )] 


(186b) 
( 186c ) 
(l86d) 


Computed values of Cf and Cf are listed in the following tabulation. 
Corresponding values of di, d 2 , Ai, and in equations (ll8) and (112) are 
also given. 
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k 


C\J 

.<VH 

lo 

dx 

d 2 

Ax 

^2 

1/6 

8.56954 

24.9568 

0.118034 

0.77900 

0.696759 

-2.12408 

1/10 

9.76336 

90.0171 

.073052 

.84874 

•790954 

- 8 . 57017 

1/20 

11.6143 

432.552 

.038646 

.97686 

. 9405 H 

-42.9388 


Results are shown in figures 4 to 8- The dimensionless magnified "boundary 
layer displacement is plotted to first and second order in figure 6 (cf. 
eq. (H4)): 


A* 



Ax + e 



(187) 


where r^y^ is the actual displacement distance corresponding to the displace- 
ment 5 = eA of the £ coordinate. These results are discussed in the next 
section. 



5 - 


0 1 L — - 1 l i j 

10 I0 2 10 3 I0 4 I0 5 I0 6 

Figure 4.- Skin friction. 
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Figure 7 .- Example of tangential velocity profile, 
k = 0.1, R = 9.251X10 3 . 



cos<£ 


Figure 8.- Example of normal velocity profile, 

k = 0 . 1 , R = 0 . 925 IXI 0 4 . 



Comparison With a Numerical Solution of the Model Problem 


The results of applying the asymptotic -expansion method to the model for 
viscous hypersonic flow described above can be compared with a numerical local 
similarity solution of the same model problem. The equations to be solved 
are (l03), (±0b) , and (105) . This problem will be seen to possess the prop- 
erty of local similarity near 9 = 0 . The solution obtained is similar to 
Probstein’s (ref. 48 ) , but for greater simplicity the viscosity has been taken 
to be constant here. 

Consider the following assumed separation of variables: 

* = ^- n2 - c P e l g (£ l (188) 

R 

(The presence of the factors l/R and e^ in equation (l88) will lead to 
elimination of those same factors from the differential equation.) For con- 
venience 9 define a parameter 


and a new variable 


N = kR(l + d) 


Pco^co r S 

LL 

*2 


5 - 


log 


1 + y 

1 + d 


(189) 


( 190 ) 


Substitution of equation (l88) into (103), (104), and (105) gives the differ- 
ential equation 


G"" - 2G m - 50" + 6G' = -2(cos cp)G(G m - 2G" - 5G* + 6G) (191) 

where ( ) * denotes differentiation with respect to either £ or £ , and the 
boundary conditions: 


At i = 0, i = £ o = -£ g = -log (1 + d): 


G = 0 

( 192 a) 

G f = 0 

(192b) 

i = 5 a s log (1 + d), f = o: 


G = (l/2)N 

( 192 c) 

G* = N(l/k - 1/2) 

(I92d) 

G" - _V f G "' - 2. + is + N (l- - 3. + 

G cos cp \ N k 2j M \k 2 k + 2 ) 

(I92e) 
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One sees that the assumed form, equation (l 88 ), did not lead to true 
separation of variables, hut that the well-known condition of local similarity 
near 9=0 is obtained by taking cos 9 equal to unity in equations (l 9 l) 
and (l92e). Then the resulting fourth-order ordinary differential equation 
can he solved with the five boundary conditions, equation ( 192 ), one of which 
is required to determine the unknown d. The two input parameters (to be 
specified) are N and k. The fourth-order differential equation ( 191 ) with 
cos 9 = 1 can be integrated once to obtain 

G m - 2G" - 5G* + 6 G = Ce S 

where 

g = - L 2G d£ = - P 2G d£ 

O 

and where 

C = (G ,, ‘ - 2G" - 5G* + 6 G)^_^ 

Since three conditions are specified at the shock and only two at the 
body, it is most convenient to start integrating at the shock, where a fourth 
condition (the value of G m or C), must be guessed and iterated. The problem 
of the unknown shock location is eliminated by using £ as the independent 
variable, which has the value zero at the shock. Equations (l92e)_ and ( 195 ) 
provide two simultaneous equations to solve for G" and G ,u at £ = 0 in 
terms of C if C is to be specified and iterated upon. 

To solve the problem numerically, one estimates a value for C, 
integrates equation ( 193 ) with decreasing £, and stops when G’ = 0. By an 

appropriate iteration scheme, one adjusts the value of C until condition 

(192a) is satisfied at the same value of £ where (l92b) is satisfied. That 
location is then denoted as £ = £ = = -log (l + d) , from which £ and d 

are determined. It is convenient to use the form ( 193 ) ^ rather than (l9l ) , 
and to iterate on values of C rather than on G ,M (^ = 0) because an initial 
estimate for C is more easily obtained. Furthermore, when N is large, C 
is very small and the iteration becomes difficult. The_diff iculty is remedied 
by estimating (from boundary-layer theory) a value of £, say £ , where the 
right side of equation (193) begins to become significant. In that case, one 
keeps the right side of equation ( 193 ) at zero until is reached, and then 

replaces the right side by qe g , where q is some specified small number and 

g = ri - 2 G dr 

Since C is taken to be zero for the first part of the integration in this 
case, one iterates instead on the value of ^ 1 , corresponding to a specified 
q, until the conditions at the body are satisfied. Different small values of 


( 193 ) 

( 194 ) 

( 195 ) 
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q are used for each case, and the results compared, to ensure that q is 
small enough that the solution obtained is not affected, to within the desired 
numerical accuracy, by this procedure. 


The dimensionless velocity components 
solution by 


(see eqs. (99)) are given in this 


u = (l/R) (sin cp)e"^(G' + G) j 
v = - 2 (l/R)(cos cp)e ^G 


(196) 


and the skin-friction coefficient, defined in equation (l86a), is 


C f = (2/kR 2 ) ( sin cp)G n (£ = 0) (197) 

The dimensionless pressure gradient along the "body surface is (with the local 
similarity condition) : 


dp *] (sin cp)G Mf (£ = 0 ) 


_c)(sin cp)j 


kR“= 


(198) 


To find the complete viscous displacement, one needs the complete outer 
solution for flow over a displaced surface at £ = 6 with no boundary layer 
present. That is, for the outer solution, the no -slip condition is replaced 
by the condition that there be no boundary layer, in line with the concept 
that £ = & is the frictionless surface that would produce the flow field 
that corresponds to the outer solution, which is approached exponentially by 
the exact solution. The complete outer solution is 4-° = (sin 2 cp)f°(£) > where 
f° is the solution e^G(£)/R that satisfies equations (191) and (192) but 
with equations (192a) and (l92b) replaced by f° = 0 at t; = S and the condi- 
tion that no rapid (exponential) variations with £ near £ = 6 are per- 
mitted. The solution for f° is found to have the same form as equation ( 125 ), 
but the constants are found from the complete boundary conditions at the shock 
and thus depend on R for v = 1 . The condition f° = 0 at £ = 5 then 
gives the following equation (refer to p. 42 ) : 


3(1 - k) 2 + 5 ( 4 k - l)(l + d*) + 2(1 - k)(l - 6k) (l + d*) ! 


, 20v 

R (1 + d) 


1 + | k)(l + d*) 2 - (l - k)(l - 6k) (l + d*) 5 


= 0 (199) 


For a given case, N is specified and d is determined, as discussed above, 
so d* is determined by this equation. Then 5 is determined from equa- 
tion (ll7) as 

5 - ^ (rrlf) (200) 

and A' = R 1/2 y & is found from equation (ll6) . 

63 



Although the calculations were made for both v = 0 and v = 1, only those 
for v = 1 are shown in figures 4 to 8. The values of Cf for v = 0 and 
v = 1 differ very little for R > 500. The differences are much less than the 
differences shown between the second-order solution and the exact numerical 
solution for v = 1. Similarly, the differences in the values of d for 
v = 0 and v = 1 were very small for R > 100. 

In the plot of Jr Cf/sin cp (fig. 4 ), first-order boundary-layer theory 
gives a constant value, shown as a dashed line for each k, Note in partic- 
ular that the second-order solution follows closely the exact numerical solu- 
tion of the model for fairly low R when k is not too small. However, 
apparently the asymptotic convergence to the exact solution of the model prob- 
lem is slowed when k becomes small. For example, for k = l/6, the first- 
order boundary-layer solution for Cf at R = 1000 is 9-8 percent low, and 
the second-order solution corrects that error to a value less than 1.6 percent 
low. For k = l/l0, the corresponding negative errors are 2 6 percent for 
first order and 4.6 percent for second order. For k = 1/20, those negative 
errors at R = 1000 are 6l percent for first order, and l6 percent for second 
order. The percent errors increase as k becomes small. This effect of 
small k slowing the asymptotic convergence can also be observed from the 
numbers given in the table following equation (l86d). As k becomes small, 
the numerical factor of the second term Cf 2 becomes large in comparison to 
that of the first term Cf^. This can be shown to be the effect of increasing 
vorticity interaction; as k becomes small, the vorticity in the external 
flow (outer solution) becomes comparable to the vorticity generated at the 
body surface. 

Similarly, the plots of shock standoff distance and boundary -layer 
displacement, and the tabulated numerical factors following equation (l86d) 
exhibit the same slowing of convergence as k becomes small. Note in partic- 
ular that the second-order boundary -layer displacement agrees well with the 
exact solution of the model problem for very low R when k = l/6. 

Figures 7 and 8 illustrate the substantial improvement of second -order 
theory in describing the u and v velocity profiles. First-order boundary- 
layer theory for the tangential velocity u goes asymptotic to a constant 
value equal to the first-order inviscid surf ace . speed. For k = 1/10 and 
R = 9.25IX10 3 , the second-order outer and inner solutions for u and v agree 
almost exactly with the exact numerical solution of the model with the inner 
solution agreeing almost exactly near y = 0, the outer agreeing more closely, 
outside the boundary layer. 


CONCLUDING REMARKS 


It has been shown how a new approach to constructing inner and outer 
expansions allows the derivation and use of a stronger form of the rule for 
matching the expansions in some singular perturbation problems of the 
boundary-layer type. By use of this "displacement matching," all displace- 
ment effects are retained explicitly, and explicit asymptotic outer boundary 
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conditions for the inner function are provided. Application of this method 
to viscous hypersonic flow over a blunt body has led to an investigation of 
the role of boundary -layer displacement, including higher -order displacements, 
not heretofore considered. In particular, second-order displacement of the 
boundary layer has been calculated, and its relationship to the shock standoff 
distance indicated. The outer and inner expansions of all functions and deri- 
vatives of the solution match according to the same stronger form of the 
matching principle. Comparison of the second -order solution, including second- 
order viscous displacement, with an exact numerical solution of the model prob- 
lem has exhibited very close agreement for Reynolds number only moderately 
large, and has exhibited a slowing of the asymptotic convergence as the shock 
density ratio k becomes small, indicating the increasing effects of 
vorticity interaction. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Nov. 10, 19 66 
124-07-02-23 


65 


I 


APPENDIX 


PRINCIPAL NOTATION 


A 


a nm 


C 

Ci 


Cf 

D* 


D 1 


D° 


d* 


-"n 


n 


^ 21^22 

f 

f* 

f 1 


see equations (154) and ( 15 6 ) 

coefficients in expansions; see equations (l 25 c) , ( 127 ), (l40b) , 
and (l42) 

constant of integration first appearing in equation ( 193 ) 

constant of integration first appearing in equation ( 162 ); 
evaluated in (l64) 

skin-friction coefficient, equations (l86) 
constants defined in equations (l86) 

domain of independent variables defined on pages 17 and 18 or on 
page 34 

domain of magnified independent variables where the magnified 
inner function is analytic with respect to (see eq. ( 60 c)) 

or where the series for F converges in equations (88) 

domain of independent variables where outer function f° is 

analytic with respect to e 0 (see eq. ( 55 ^)) or where the series 
for f° converges in equations ( 87 ) 

r s - rt 

dimensionless shock standoff distance, 

r b 

defined on page 42, determined by equation (l99) 
see equations (108), (llO), and (ll8) 
terms of expansion of d, equation (llO) 
magnified dependent variable, equation ( 60 c) 

functions of x,Y in equations ( 62 ) and ( 92 b); in the blunt -body 
problem, functions of tj; see equations ( 152 ) and ( 159 ) 

defined by equations ( 165 ) through ( 167 ) 

generally a function of x,y, z,e (see pages 5 and 17 ) 

defined on pages 17 , 18; also see page 34 

inner function (eqs. ( 60 b,c); also eqs. (88)) 
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f° outer* function (see pp. 21, 22, and 3 4) 

f n generally functions of x,y in an asymptotic expansion (see p. 17) 

f° outer displacement dependent variable ; see equations (68), (69), and 
(90b) 

f n functions of x,y in the terms of an expansion, equations (69), (91); 
in the blunt body problem, see equations (123), (138) 

G(£) see equation (l88) 

g arbitrary function; and function of £ defined in equation (19^) tor 

the blunt body problem 

h specific enthalpy 

K. constants defined on page 53 

J 

k Po > 2 

L linear differential operator defined by equation (163) 

N Reynolds number defined in equation (189) 

P magnified inner pressure function defined on page 55 

P n terms of the expansion (174) 
p pressure 

p dimensionless pressure function (p. 38) 

pi inner pressure function defined on page 5^ 

— o 

p outer pressure function 

q* defined on page 18 

q 1 defined in equations (60b), (60c), (88); see also equations (132) and 

(13*0 

q° defined on page 21; also equations (107) and (ill) 

R Reynolds number, p^V^r^/q^, equation (lOl) 

r radial distance measured from sphere center 

r§ defined on page 42 

u,v dimensionless velocity components; see page 38 
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flow velocity vector and magnitude of velocity, respectively 
defined on page 22; see also pages 25 and 2 6 

coordinate x, used along with y in the transformation to displace- 
ment variables (eqs. ( 67 )) 

magnified independent variable, equation ( 60 a) 

general independent variable; in the blunt body problem y is 
dimensionless and r^y is the distance from the body surface 

e 5 - 1 (eq. (ll4)); r-^yg is the displacement distance from the body 
surface (see also p. 58 ) 

independent displacement variable, y - S; see pages 27, 28 , and 35 

magnified independent variable in the blunt body problem (eq. ( 13 I)) 

parameter in example problems 2 (p. 7 ) and 3 (p. 9) > also used on 
page 31 

see equations ( 157 ), ( 158 )., (l6l) 
magnified displacement, 8/crg 

terms of expansion of A, equations (66) and (112) 

dimensionless magnified displacement of y coordinate, equation (ll4) 

general dimensionless displacement in boundary-layer type problems; 
see equations (22c), (24b), (34b), (591), (64), and (112) 

generally a small parameter (see p. 5)j in the blunt body problem, 
e = R -1/2 (eq. (137)) 

see page 25 

see page 21 

independent variable in blunt body problem (eqs. (lOO)) 
displacement variable, £ - B 
defined in equation ( 190 ) 

independent variable for boundary -layer solution; see equation ( 1531 ) 
see equations (l53c), (156) 



[X parameter in example problem 3 (p. 8); also viscosity coefficient 

in blunt body problem 

v artificially inserted parameter, first appearing in condition (l 05 c); 

discussed on pages 40, 41 

(■ independent variable for the boundary -layer solution; see 

equations (l53a) and (156) 

p mass density 

a i order of magnitude of width of nonuniform region (see eq. ( 60 a) and 

pp. 46, 47) 

0*2 order of magnitude of dependent variable in region of nonuniformity 

(see eq. ( 60 c) and pp. 46, 47) 

order of magnitude of 5; usually = a± (see eq. ( 65 )) 

o(cp) function of integration defined on page 50 

cp angle measured from stagnation streamline, centered at sphere center 

(fig. 3) 

cp variable cp, used along with £ in transformation to displacement 

variables (eq. (ll9a)) 

1 magnified inner stream function in the blunt body problem, 

equation ( 131 ) 

terms of expansion ( 133 ) 

\|r dimensionless Stokes stream function defined by equations (99) 

\|r n terms of expansion ( 109 ) 

i]/ 1 2 defined on page 46 

\| r° defined on page 41 

-ty° defined in equation (ll9b) 

\jr n terms of expansion of ijr 0 , equation (ll9c) 

w dimensionless vorticity, see page 45 

Subscripts 

2 value behind shock on stagnation streamline 

co value in free stream 


69 



b value at body surface or boundary 
n index denoting term of an expansion 

s value behind shock 

cp, £ denote partial differentiation with respect to cp or respectively 

cp,Z denote partial differentiation with respect to cp or Z, respectively 

Superscripts 

i inner function or function transformed from outer displacement 

variables to inner magnified variables 

o outer function or function transformed to outer displacement 
variables 
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